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1. Introduction and statement of the results 

In this paper we show that for a class of open quantum systems satisfying a natural 
dynamical assumption (see §2.21) the study of the resolvent, and hence of scattering, and of 
resonances, can be reduced, in the semiclassical limit, to the study of open quantum maps, 
that is of finite dimensional quantizations of canonical relations obtained by truncation of 
symplectomorphisms derived from the classical Hamiltonian flow (Poincare return maps). 

We first explain the result in a simplified setting. For that consider the Schrodinger 
operator 



P{h) = -h^A + V{x) - 1 



V e 



and let be the corresponding classical flow on T*]R" 3 {^lO- 



x'{t) = 2e(t) 



$*(x,o = (x(t),eW), 

e{t) = -dVixit)), x(0) 



X . 




Figure 1. An example of a potential, V G C^(M^), to which the results 
apply: the Hamiltonian flow is hyperbolic on the trapped set in a range of 
energies - see [381 Appendix c]. In this example each energy surface p~^{E) 
is three dimensional, so the Poincare section is two dimensional as shown in 
Fig.H 
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Equivalently, this flow is generated by the Hamilton vector field 




associated with the classical Hamiltonian 



(1.3) p(a;,e) = |eP + na^)-l. 



The energy shift by —1 allows us to focus on the quantum and classical dynamics near the 
energy E' = 0, which will make our notations easieiQ. We assume that the Hamiltonian 
fiow has no fixed point at this energy: dp\p-i(Q)^ 0. 

The trapped set at any energy E is defined as 

(1.4) Ke = {(x, e r*R" : p(x, = E, $*(x, remains bounded for all t G M} . 

The information about spectral and scattering properties of P = P{h) in fll.ip can be 
obtained by analyzing the resolvent of P, 



and its meromorphic continuation - see for instance [33] and references given there. More 
recently semiclassical properties of the resolvent have been used to obtain local smoothing 
and Strichartz estimates, leading to applications to nonlinear evolution equations - see [H] 
for a recent result and for pointers to the literature. In the physics literature the Schwartz 
kernel of R{z) is referred to as Green's function of the potential V . 

The operator P has absolutely continuous spectrum on the interval [—1, oo); nevertheless, 
its resolvent R{z) continues meromorphically from Im z > to the disk -D(0, 1), in the sense 
that xP(z)x, X ^ C^(]R"), is a meromorphic family of operators, with poles independent 
of the choice of x ^ (see for instance [HI Section 3] and [39], Section 5]). 

The multiplicity of the pole z & D{Q, 1) is given by 



where the integral runs over a sufficiently small circle around z. 

We now assume that at energy = 0, the fiow is hyperbolic on the trapped set Kq 
and that this set is topologically one dimensional. Hyperbolicity means [2H Def. 17.4.1] 
that at any point p = (x, ^) G -f^o the tangent space to the energy surface splits into the 



R{z) = {P-z) 



-1 



Im z > , 





""^There is no loss of generality in this choice: the dynamics of the Hamiltonian + V{x) at some energy 
i? > is equivalent with that of + V /E —1 at energy 0, up to a time reparametrization by a factor \/E. 
The same rescaling holds at the quantum level. 
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Figure 2. A schematic view of a Poincare section S = UjSj for Ke inside 
p~^{E). The flow near Ke can be described by an ensemble of symplectomor- 
phisms between different components - see §2.21 for abstract assumptions 
and a discussion why they are satisfied when the flow is hyperbolic on Ke 
and Ke has topological dimension one. The latter condition simply means 
that the intersections of Ke with E/s are totally disconnected. 

this decomposition is preserved through the flow, and is characterized by the following 
properties: 

(1.6) 3C>0, 3A>0, \\dexptHp{p)v\\ <Ce-^^^^\\v\\, \/vEEj,±t>0. 

When Kq is topologically one dimensional we can find a Poincare section which reduces the 
flow near Kq to a combination of symplectic transformations, called the Poincare map F: 
see Fig J2] for a schematic illustration and §2.21 for a precise mathematical formulation. The 
structural stability of hyperbolic flows Thm. 18.2.3] implies that the above properties 
will also hold for any energy ii^ in a sufficientlys short interval [—6, 6] around E = 0, in 
particular the flow near Ke can be described through a Poincare map Fe- 

Under these assumptions, we are interested in semiclassically locating the resonances of 
the operator P{h) in a neighbourhood of this energy interval: 

n{S,Mo,h) = [-6,6]+t[-Mohlogil/h),Mohlog{l/h)] , 

where S, Mq are independent of /i G (0, 1]. Here the /ilog(l//i)-size neighbourhood is natural 
in view of results on resonance free regions in case of no trapping - see |26] . 

To characterize the resonances in 7l{6, Mq, h) we introduce a family of "quantum prop- 
agators" quantizing the Poincare maps Fe- 
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Theorem 1. Suppose that $* is hyperbolic on Kq and that Kq is topologically one dimen- 
sional. More generally, suppose that P{h) and $* satisfy the assumptions of ^2.1\ - ^2.2 . 

Then, for any 6 > small enough and any Mq > 0, there exists ho > such that there 
exists a family of matrices, 

{M{z,h), z en{6,Mo,h), he{0,ho]}, 

holomorphic in the variable z, and satisfying 

/i-'^+VCo < rank M{z, h) < Coh-''+^ , Co > 1 , 

such that for any h G (0, Hq], the zeros of 

({z, h) = det(J-M(2,/i)), 

give the resonances of P{h) in 7l{6, Mq, h) , with correct multiplicities. 

The matrices M{z, h) are open quantum maps associated with the Poincare maps Fn^z de- 
scribed above: for any L > 0, there exist a family of h-Fourier integral operators, {Ai{z, h)}, 
quantizing the Poincare maps F^cz (see ^2.3.2\ and ^3. 3\) . and projections Uh (see ^5.2.2\) 
of ranks 

/i-'^+VCo < rankn^, < Coh~''+^ , 

such that 

(1.7) M{z,h) =U,,M{z,h)Uh + 0{h^) . 

The statement about the multiphcities in the theorem says that 

, , 1 fC{w), 



— tr /(/ - M(w))-^M'(w)dw . 
2vr2 Jz 



A more precise version of Theorem [H involving complex scaling and microlocally deformed 
spaces (see §3.41 and §3.51 respectively), will be given in Theorem [2] in §5.41 In particular 
Theorem [2] gives us a full control over both the cutoff resolvent of P, xR{z)Xj ^ind the 
full resolvent (-Pg,/? — of the complex scaled operator Pq^r, in terms of the family of 
matrices M{z, h); for this reason, the latter is often called an effective Hamiltonian for P. 

The mathematical applications of Theorem[T]and its refined version below include simpler 
proofs of fractal Weyl laws |13] and of the existence of resonance free strips [HT]. The 
advantage lies in eliminating flows and reducing the dynamical analysis to that of maps. 
That provides an implicit second microlocalization without any technical complication (see 
[13| §5]). The key is a detailed understanding of the operators A4{z,h) stated in the 
theorem. 



Relation to semiclassical trace formulcB. The notation C{z., h) in the above theorem hints at 
the resemblance between this determinant and a semiclassical zeta function. Various such 
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functions have been introduced in the physics hterature, to provide approximate ways of 
computing eigenvalues and resonances of quantum chaotic systems - see |171 |20l [10] . 

These semiclassical zeta functions are defined through formal manipulations starting 
from the Gutzwiller trace formula - see |l2] for a mathematical treatment and references. 
They are given by sums, or Euler products, over periodic orbits where each term, or factor 
is an asymptotic series in powers of h. Most studies have concentrated on the zeta function 
defined by the principal term, without /i-corrections, which strongly resembles the Selberg 
zeta function defined for surfaces of constant negative curvature. However, unlike the case 
of the Selberg zeta function, there is no known rigorous connection between the zeroes 
of the semiclassical zeta function and the exact eigenvalues or resonances of the quantum 
system, even in the semiclassical limit. Nevertheless, numerical studies have indicated that 
the semiclassical zeta function admits a strip of holomorphy beyond the axis of absolute 
convergence, and that its zeroes there are close to actual resonances [TU| 

The traces of M{z,h)'', A; G N admit semiclassical expressions as sums over periodic 
points, which leads to a formal representation of 



as a product over periodic points. That gives it the same form as the semiclassical zeta 
functions in the physics literature. In this sense, the function ({z, h) is a resummation of 
these formal expressions. As will become clear from its construction below, the operator 
M{z, h) is not unique: it depends on many choices which affect the remainder term 0{h^) 
in fll.7l) . However, the zeroes of ({z,h) in Tl{S, Mq, h) are the exact resonances of the 
quantum Hamiltonian. 

Comments on quantum maps in the physics literature. Similar methods of analysis have 
been introduced in the theoretical physics literature devoted to quantum chaos. The clas- 
sical case involves a reduction to the boundary for obstacle problems: when the obstacle 
consists of several strictly convex bodies, none of which intersects a convex hull of any other 
two bodies, the fiow on the trapped set is hyperbolic. The reduction can then be made 
to boundaries of the convex bodies, resulting with operators quantization Poincare maps 
- see Gaspard and Rice [12] , and for a mathematical treatment Gerard [TB] , in the case of 
two convex bodies, and [2HI §5.1], for the general case. FigO illustrates the trapped set in 
the case of three discs. The semiclassical analogue of the two convex obstacle, a system 
with one closed hyperbolic orbit, was treated by Gerard and the second author in [19]. The 
approach of that paper was also based on the quantization of the Poincare map near this 
orbit. 

A reduction of a more complicated quantum system to a quantized Poincare map was 
proposed in the physics literature. Bogomolny [1] studied a Schrodinger operator P{h) 
with discrete spectrum, and constructed a family of energy dependent quantum transfer 
operators T{E, h), which are integral operators acting on a hypersurface in the configuration 




k=l 
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Figure 3. This figure, taken from [31], sfiows tlie case of symmetric three 
disc scattering problem (left), and the associated Poincare section (right). 
The section is the union of the three coball bundles of circle arcs (in red) 
parametrized by s (the length parameter on the circle, horizontal axis), and 
cos if (vertical axis), where if is the angle between the velocity after impact 
and the tangent to the circle. Green, blue,red strips correspond to different 
regions of forward escape; they are bounded by components of the stable 
manifold. The trapped set, T, shown in yellow, is the intersection of the 
latter with the unstable manifold. 

space. These transfer operators are asymptotically unitary as /i — )■ 0. The eigenvalues of 
P{h) are then obtained, in the semiclassical limit, as the roots of the equation det(l — 
T{E)) = 0. Smilansky and co-workers derived a similar equation in the case of closed 
Euclidean 2-dimensional billiards [13], replacing T{E) by a (unitary) scattering matrix 
S{E) associated with the dual scattering problem. Prosen [53] generalized Bogomolny's 
approach to a nonsemiclassical setting. Bogomolny's method was also extended to study 
quantum scattering situations [T71 152] . 

Open quantum maps have first been defined in the quantum chaos literature as toy 
models for open quantized chaotic systems (see [291 §2-2], [30], §4.3] and references given 
there). They generalized the unitary quantum maps used to mimic bound chaotic systems 
[TT] . Some examples of open quantum maps on the 2-dimensional torus or the cylinder, 
have been used as models in various physical settings: Chirikov's quantum standard map 
(or quantum kicked rotator) was first defined in the context of plasma physics, but then 
used as well to study ionization of atoms or molecules [H], as well as transport properties 
in mesoscopic quantum dots [IHj. Other maps, like the open baker's map, were introduced 
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as clean model systems, for which the classical dynamics is well understood [361 EHj • The 
popularity of quantum maps mostly stems from the much simplified numerical study they 
offer, both at the quantum and classical levels, compared with the case of Hamiltonian flows 
or the corresponding Schrodinger operators. For instance, the distribution of resonances 
and resonant modes has proven to be much easier to study numerically for open quantum 
maps, than for realistic flows [3 (33 1221 [23 127] . Precise mathematical definitions of quantum 
maps on the torus phase space are given in [291 §4.3-4.5]. 



Organization of the paper. In the remainder of this section we give assumptions on the 
operator P and on the corresponding classical dynamical system, in particular we introduce 
a Poincare section "E and map associated with the classical flow. We refer to results 
of Bowen and Walters [8] to show that these assumptions are satisfied if the trapped 
set supports a hyperbolic flow, and is topologically one dimensional, which is the case 
considered in Theorem [H 

In ^we recall various tools needed in our proof: pseudodifferential calculus, the concept 
of semiclassical microlocalization, local /i- Fourier integral operators associated to canonical 
tranformation (these appear in Theorem [T|), complex scaling (used to define resonances 
as eigenvalues of nonself adjoint Fredholm operators), microlocally deformed spaces, and 
Grushin problems used to define the effective Hamiltonians. 

In §1] we follow a modified strategy of [12] and construct a microlocal Grushin problem 
associated with the Poincare map on S. No knowledge of that paper is a prerequisite but 
the self-contained discussion of the problem for the explicit case of given in [121 §2] can 
illuminate the complicated procedure presented here. In [121 §2] one finds the proof of the 
classical Poisson formula using a Grushin problem approach used here. 

Because of the hyperbolic nature of the flow the microlocal Grushin problem cannot 
directly be made into a globally well-posed problem - see the remark at the end of §11 This 
serious difficulty is overcome in §S] by adding microlocal weights adapted to the flow. This 
and suitably chosen finite dimensional projections lead to a well posed Grushin problem, 
with an effective Hamiltonian essentially given by a quantization of the Poincare map: this 
fact is summarized in Theorem [21 from which Theorem [1] is a simple corollary. 
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2. Assumptions on the operator and on classical dynamics 
Here we carefully state the needed assumptions on quantum and classical levels. 

2.1. Assumptions on the quantum Hamiltonian P{h). Our results apply to operators 
P{h) satisfying general assumptions given in pTJ §3.2] and [l3l (1.5), (1.6)]. In particular, 
they apply to certain elliptic differential operators on manifolds X of the form 

J 

X = X^U|J (M"\5Mn(0,i?)) , 
i=i 

where i? > is large and is a compact subset of X. The reader interested in this higher 
generality should consult those papers. 

Here we will recall these assumptions only in the (physical) case of differential operators 
on X = M". We assume that 

(2.1) p{h) = j2(^c.{x,h){hD,r, 

\a\<2 

where a„(x, h) are bounded in C°°{W), a^ix, h) = a° (x) + C(/i) in C°°, and 

is independent of /i for |a| = 2. Furthermore, for some Cq > the functions a^ix, h) have 

holomorphic extensions to 

(2.2) {xgC" : |Rea;|>Co, | Imx| < | Rea;|/Co} , 

they are bounded uniformly with respect to h, and aa{x, h) = a° (x) + 0{h) on that set. 

Let P{x,^) denote the (full) Weyl symbol of the operator P, so that P = P'^{x] hD; h), 
and assume 

(2.3) p(a;,e;/,)->e2_^ 

when X — oo in the set fl2.2p . uniformly with respect to {^,h) G A'x]0, 1] for any compact 
set (E (here, and below, d means that the set on the left is a pre-compact subset of 
the set on the right). We also assume that P is classically elliptic: 

(2.4) p2{x, = Yl «"(^)^" ^ \ {0}' 

\a\=2 

and that P is self-adjoint on //^(M") with domain H'^{W^). The Schrodinger operator (11. ip 
corresponds to the choices ^|ct|=2'^a^" ~ I^P' = for |a| = 1, and ao(a^) = V{x) — 1. 
The assumption (12. 3 p show that we can also consider a slowly decaying potential, as long 
as it admits a holomorphic extension in (12. 2p . 
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2.2. Dynamical Assumptions. The dynamical assumptions we need roughly mean that 
the flow $* on the energy shell p~^(0) C T*X can be encoded by a Poincare section, 
the boundary of which does not intersect the trapped set Kq. The abstract assumptions 
below are satisfled when the flow is hyperbolic on the trapped set which is assumed to be 
topologically one dimensional - see Proposition 12.11 

To state the assumption precisely, we notice that 

(2.5) P(x,0= 

|q|<2 

is the semi-classical principal symbol of the operator P{x, hD; h). We assume that the 
characteristic set of p (that is, the energy surface p~^(0)) is a simple hypersurface: 

(2.6) dp y^O on p~\0). 
Like in the introduction, we denote by 

$* = exp{tHp) : T*X T*X 

the flow generated by the Hamilton vector fleld Hp (see (11. 2p ). 

Our assumptions on p{x,^) ensure that, for E close to 0, we still have no flxed point in 
p~^{E), and the trapped set Ke (deflned in (II. 4p ) is a compact subset of p~'^{E). 

We now assume that there exists a "nice" Poincare section for the flow near Kq, namely 
flnitely many compact contractible smooth hypersurfaces Sfc C p~^{0), k = 1,2, . . . , N with 
smooth boundaries, such that 

(2.7) 9Sfcni^o = 0, SfcnSfc. = 0, k^k', 

(2.8) Hp is transversal to uniformly up to the boundary. 

For every p E Kq, there exist p_ G ^j_{p) , p+ G ^j+{p) 

(2.9) of the form p± = $^*±('')(p), with < t±(p) < W < oo, such that 

{$*(p); -t_(p) <t <t+(p), t ^OjnSfc = 0, \/k. 
We call Poincare section the disjoint union 

The functions p ^ p±{p) , p ^ t±{p) are uniquely deflned {p±{p) will be called respectively 
the successor and predecessor of p). They remain well-deflned for p in some neighbourhood 
of Kq in p~^{0)) and, in such a neighbourhood, depend smoothly on p away from S. In 
order to simplify the presentation we also assume the successor of a point p G belongs 
to a different component: 

(2.10) If p G Sfc n Kq for some k, then p+(p) G fl i^o for some I ^ k. 
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The section can always be enlarged to guarantee that this condition is satisfied. For in- 
stance, for Kq consisting of one closed orbit we only need one transversal component to 
have fl2.7p -( l2l8l) : to fulfill fl2.10p a second component has to be added. 

We recall that hypersurfaces in p~^{0) that are transversal to Hp are symplectic. In fact, 
a local application of Darboux's theorem (see for instance [231 §21.1]) shows that we can 
make a symplectic change of variables in which p = C,n and Hp = dx„- If S C {^n = 0} 
is transversal to d^^, then (xi, ■ ■ -Xn-ijC,!, ■ ■ ■ ,^n-i) can be chosen as coordinates on S. 
Since u\p-i(^Q)= X]j=i d-^j ^ ^^at means that u\j] is nondegenerate. The local normal 
form p = will be used further in the paper (in its quantized form). 

The final assumption guarantees the absence of topological or symplectic peculiarities: 
There exists a set d T*M"~^ with smooth boundary, and a symplectic 
(2.11) diffeomorphism : — t- which is smooth up the boundary together 

with its inverse. We assume that extends to a neighbourhood of Sfc in T*M.'^. 

In other words, there exist symplectic coordinate charts on S^, taking values in E^. 

The following result, due to Bowen and Walters [8], shows that our assumptions are 
realized in the case of 1-dimensional hyperbolic trapped sets. 

Proposition 2.1. Suppose that the assumptions of ^2.1\ hold, and that the flow 

is hyperbolic in the standard sense of Eqs. M.5^1.0\) . Then the existence of S satisfying 

fl2.7p - fl2.lip is equivalent with Kq being topologically one dimensional. 

Remark. Bowen and Walters [8] show more, namely the fact that the sets {S^} can be 
chosen of small diameter, and constructed such that S fl Kq forms a Markov partition for 
the Poincare map. Small diameters ensures that (12. lip holds, while, as mentioned before, 
fl2.10p can always be realized by adding some more components. 

Proposition 12.11 shows that the assumptions of Theorem [T] imply the dynamical assump- 
tions made in this section. The proof of Appendix c] shows that the following example 
of "three-bumps potential" , 

3 

P = -h^A + V{x) -1, a;eR^ V{x) = 2^exp{-R{x - Xkf) , 

k=l 

Xk = (cos(27r/i;/3), sin(27rA;/3)) , 
satisfies our assumptions as long as i? > 1 is large enough (see Fig. [1]). 

2.3. The Poincare map. Here we will analyze the Poincare map associated with the 
Poincare section discussed in §2.21 and its semiclassical quantization. 
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Figure 4. Schematic representation of the components Fik of the Poincare 
map between the sets -Dj^ and Ai^ (horizontal/vertical elhpses). The reduced 
trapped set 71 is represented by the black squares. The unstable/stable 
directions of the map are the horizontal/vertical dashed hues. 



2.3.1. Classical analysis. The assumptions in §2.21 imply the existence of a symplectic re- 
lation, the so-called Poincare map on S. 

More precisely, let us identify S^'s with using given in (12. lip , so that the Poincare 
section 



N N N 



Let us call 



The map 



k=l k=l k=l 

T == -R'o n S = I 1 7fc the reduced trapped set. 



(see the notation of (12. 9p ) is the Poincare map for ^*\ko- It is a Lipschitz bijection. The 
decomposition T = Tk allows us to define the arrival and departure subsets of T: 

Afc = {p G 71 C : p_(p) G = ?- n /(Tfc) = f{V,k) , 

For each k we call J+{k) C {1, . . . , N} the set of indices i such that is not empty (that 
is, for which 71 is a successor of Tk)- Conversely, the set J-{i) refers to the predecessors of 
%. 

Using this notation, the map / obviously decomposes into a family of Lipschitz bijections 
fik '■ T>ik — )■ Aik- Similarly to the maps p±, each fik can be extended to a neighbourhood 
of "Djfc, to form a family of local smooth symplectomorphisms 

Fik '■ — > Fik{Dik) = Aik, 
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Figure 5. Trajectories linking the boundaries of the departure set Dn. C 
and the arrival set Aik C Ej. Note the stretching and contraction implied 
by hyperbolicity. These trajectories and Dik U Aik form the boundary of the 
tube Tik defined by ([212]). 

where Dik (resp. Aik) is a neighbourhood of Vik in Sfc (resp. a neighbourhood of Aik in Sj). 
Since our assumption on Kq is equivalent with the fact that the reduced trapped set T is 
totally disconnected, we may assume that the sets {Dik}i^j_^_(k) (resp. the sets {^jfc}fceJ-(j)) 
are mutually disjoint. We will call 

Dk = ^i<^j+{k)Dik, Ai = L\k<^j_(k)Aik ■ 

Notice that, for any index i, the sets Di, Ai both contain the set 71, so they are not disjoint. 
We will also define the tubes Tik C T*X containing the trajectories between Dik and Aik'. 

(2.12) Tik = {$*(p), : P e Dik, 0<t< t+{p)} . 

See Fig. H] for a sketch of these definitions, and Fig. [5] for an artistic view of Tik 

The maps Fik will be grouped into the symplectic bijection F between Dk and Ak- 
We will also call F the Poincare map, which can be viewed as a symplectic relation on S. 
We will sometimes identify the map Fik with its action on subsets of T*]R"~^. 

Fik = K,~'^ o Fik o Kk : Dik — >Ak, Dik'^= K^'^iDik) , Aik = k^^ (Aik) . 

Using the continuity of the flow $*, we will show in §4.1.11 that the above structure can 
be continuously extended to a small energy interval z G [—6, 6] . The Poincare map for the 
flow in p~^{z) will be denoted by F^ = {Fik^z)i<i,k<N (see §4.1.1l for details). 
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In the case of Ko supporting a hyperbolic flow, a structural stability of holds in a 
stronger sense: the flows $*titz cind ^*\ko are actually orbit-conjugate (that is, conjugate up 
to time reparametrization) by a homeomorphism close to the identity. |2H Thm. 18.2.3]. 

2.3.2. Quantization of the Poincare map. In this section we make more explicit the operator 
Ai{z, h) used in Theorem [H The semiclassical tools we are using will be recalled in ^ 

Let us flrst focus on a single component Fik : -Djfc — ^ ^ifc of the Poincare map. A quantiza- 
tion of the symplectomorphism Fi^ (more precisely, of its puUback Fj^) is a semiclassical (or 
h-) Fourier integral operator, that is a family of operators M.ik{h) : L^(M"~^) — > L^(M""^), 
h G (0, 1], whose semiclassical wavefront set satisfles 

(2.13) WF'JMfc) <^ Aa: X Afc , 

and which is associated with the symplectomorphism Fik. {h-FlOs are deflned in §3.3[ and 
WF';, is deflned in below). 

Being associated to the symplectic map F^k means the following thing: for any a G 
C^{Aik), the quantum operator Op]^(a) transforms as follows when conjugated by M.ik{h) 

(2.14) M.k{hrOvl{a)Mik{h) = OpH^^.^la) + h'-'' Op]^(6) , 

where the symbol aik G ^^(T*]^'""^) is independent of a, Oik = 1 on some neighbourhood 
of Tfc in Sfc, and b G Ss{T*W^~^), for every 6 > 0. Here Op^ denotes the semiclassical Weyl 
quantization on (see eq. fl3.ip ). and Ss{T*'K"'~^) is the symbol class deflned in §3.11 

The necessity to have (5 > in fl2.14p comes from the slightly exotic nature of our Fourier 
integral operator, due to the presence of some mild exponential weights - see §3.51 below. 

The property fl2.14p . which is a form of Egorov's theorem, characterizes Aiik{h) as a 
semiclassical Fourier integral operator associated with Fik (see [121 Lemma 2] and [151 
§10.2] for that characterization). 

We can then group together the Aiik{h) into a single operator- valued matrix (setting 
Mik{h) = when i ^ J+{k)): 

Mih) : L'iR-~Y ^ L'iR-^'', Mih) = (A^.fc(/^))i<,,,<^ • 

We call this A4.{h) a quantization of the Poincare map F. 

The operators (2;, h) in Theorem [1] will also holomorphically depend on 2; G 'R-{5, Mq, h), 
such that for each z G 7^(5, Mq, /i) fl M the family {J^{z, h))he{o,i] is an /i-Fourier integral 
operator of the above sense. 

Comment on notation. Most of the estimates in this paper include error terms of the 
type 0{h°°), which is natural in all microlocal statements. To simplify the notation we 
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adopt the following convention (except in places where it could lead to confusion): 

u = v \\u - v\\ = 0(h°^)\\u\\ , 

. - II II V /II in 

< ||Tn|| + ||t;|| ^ <0(l)(||Tn|| + ||i;||) + 0(/i°^)||n||, 

with norms appropriate to context. Since most estimates involve functions u micro localized 
to compact sets, in the sense that, u — x{x,hD)u G /i°°=5^(M"'), for some x ^ C^{T*R"'), 
the norms are almost exclusively norms, possibly with microlocal weights described in 

mM 

The notation u = Ov{f) means that \\u\\y = 0{f), and the notation T = Oy^wif) 
means that ||Tu||vk = Also, the notation 

neigh(yl, B) for Ac B, 

means an open neighbourhood of the set A inside the set B. 

Starting with §U we denote the Weyl quantization of a symbol a by the same letter 
a = a^{x, hD). This makes the notation less cumbersome and should be clear from the 
context. 

Finally, we warn the reader that from §1] onwards the original operator P is replaced by 
the complex scaled operator Po^r, whose construction is recalled in §3.41 Because of the 
formula f l3.13p . that does not affect the results formulated in this section. 



3. Preliminaries 

In this section we present background material and references needed for the proofs of 
the theorems. 

3.1. Semiclassical pseudodifferential calculus. We start by defining a rather general 
class of symbols (that is, /i-dependent functions) on the phase space T*M'^. For any S G 
[0,1/2] and m,k eR, let 

S^^\T*R'^) = {ae C°°(T*M'' X (0,1]) : Va G N^ /3 G N^ 3C„^ > 0, 

where (0 = (1 + 1^1')^- 

Most of the time we will use the class with 5 = in which case we drop the subscript. 
When m = A; = 0, we simply write S{T*R'^) or S for the class of symbols. In the paper 
d = n (the dimension of the physical space) ot d = n — 1 (half the dimension of the Poincare 
section), and occasionally (as in 02.131) ) d = 2n — 2, depending on the context. 
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The quantization map, in its different notational guises, is defined as follows 

a^u = Op^{a)u{x) = a"'(x, hD)u{x) 

and we refer to [121 Chapter 7] for a detailed discussion of semiclassical quantization (see 
also [ini Appendix]), and to [T51 Appendix D.2] for the semiclassical calculus for the symbol 
classes given above. 

We denote by \Ef™''^(M'^) or '^"^'^(M.'^) the corresponding classes of pseudodifferential op- 
erators. The quantization formula fl3.ip is bijective: each operator A G \E'™''^(]R'^) is exactly 
represented by a unique (full) symbol a{x,^; h). It is useful to consider only certain equiv- 
alence classes of this full symbol, thus defining a principal symbol map - see [151 Chapter 



The combination o Op^ is the natural projection from S^'^ onto 5'^>^y'5'™-i'^-i+2'5_ 
main property of this principal symbol map is to "restore commutativity" : 

ah{A o B) = ah{A)(ThiB) . 

Certain symbols in S'™'°(T*]R'^) admit an asymptotic expansion in powers of h, 

(3.2) a{x,^;h) ~ ''^^h^ aj{x,^), aj E 5"""-''° independent of h, 

such symbols (or the corresponding operator) are called classical, and make up the subclass 
5'™'*^(T*M'^) (the corresponding operator class is denoted by \E'^'°(]R'^)). For any operator 
A E \Ef™'°(]R'^), its principal symbol ^^(A) admits as representative the /i- independent 
function ao(x, ^), first term in (13.21) . The latter is also usually called the principal symbol 
of a. 

In §3.51 we will introduce a slightly different notion of leading symbol, adapted to a 
subclass of symbols in S'(T*M) larger than Sci{T*'R'^). 

The semiclassical Sobolev spaces, if^(M'^) are defined using the semiclassical Fourier 
transform, 

(3.3) lkll?,;= / {O'-lJ'MOl'd^ , -^^"(^) = T^TTW^ / ^i^)e-'^''^^^'dx . 

Unless otherwise stated all norms in this paper, || • ||, are norms. 

We recall that the operators in \E'(]R'') are bounded on uniformly in h, and that they 
can be characterized using commutators by Beals's Lemma (see [121 Chapter 8] and [^ 
Lemma 3.5] for the Ss case): 

(3.4) Ae^,(x) ^ I "/'^^-'t-f "-fr f^"'"":r' . 

^ ^ I for Imear functions «j(x, 4) on M'^ X M", 
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where ad^A = [5, A]. 

For a given symbol a G S{T*W^) we follow [12] and say that the essential support is 
contained in a given compact set K d T*M'^, 

ess-supp;, aC K T*R'^ , 

if and only if 

Vx G S{T*R'^) , suppx n = =^ xae /i~^(T*M^) . 
The essential support is then the intersection of all such K^s. 

Here ^ denotes the Schwartz space. For A G \I^(]R'^), A = Op'^{a), we call 
(3.5) WFh{A) = ess-supp;j a. 

the semiclassical wavefront set of A. (In this paper we are concerned with a purely semi- 
classical theory and will only need to deal with compact subsets of T*M'^. Hence, we won't 
need to define noncompact essential supports). 

3.2. Microlocalization. We will also consider spaces of functions (strictly speaking, 
of /i-dependent families of functions) which are microlocally concentrated in an open set 
V d T*R'^: 

H{V) = {u = {u{h) G L2(M'^))^g(o i], such that 

(3-6) 3a>0, V/iG (0,1], ||n(/i)|U2(K.) <a, 

3 X G C^{V) , x'^ix, hD,) u{h) = u{h) + Oy{h°^)] . 

The semiclassical wave front set of n G H{V) is defined as: 
(3.7) 

WFh{u) =T*R'^\{{x,^) eT*R'^ : 3aeS{T*R'^), a(x,0 = l, \\a'" u\\l2 = 0{h'^)} . 

The condition (13.71) can be equivalently replaced with u = Oy{h°°), since we may always 
take a G y(T*R'^). This set obviously satisfies WFh{u) d V. Notice that the condition 
does not characterize the individual functions u{h), but the full sequence as /i — )■ 0. 

We will say that an /i-dependent family of operators T = {T{h))he(o,i] ■ ^(K"') ^ y"{R'') 
is semiclassically tempered if there exists L > such that 

||(a;)"^T(/i)u||^-L < C/i"^||(x)^n||^L, /i G (0, 1) , {x) = {I + x^f'^ . 

Such a family of operators is microlocally defined on V if one only specifies (or consid- 
ers) its action on states u G H{y), modulo 0,yi^,y{h'^). For instance, T is said to be 
asymptotically uniformly bounded on H{y) if 

(3.8) 3Cr>0VnGi/(\/) >0, V/iG (0,/iT,«), \\T{h)u{h)\\L2(^^u^ < CtC^ . 

Two tempered operators T, T' are said to be microlocally equivalent on V , iff for any 
u G H{V) they satisfy ||(T - T')n||i2(Rfc) = 0(/i~); equivalently, for any x e C,°°(V^), 
\\{T-T')xnL^^L- = 0{h^). 



FROM OPEN QUANTUM SYSTEMS TO OPEN QUANTUM MAPS 17 

If there exists an open subset W (s T*M^ and L G M such that T maps any u G H(y) 
into a state Tu G H{W), then we will write 

T = T{h) : ^ , 

and we say that T is defined microlocally in W x V. 

For such operators, we may define only the part of the (twisted) wavefront set which is 
inside W xV: 
(3.9) 

WF';,(r) n{W xV) = {W xV)\ {(p',p) eW xV ■. 3a G S{T*R^), b G S{T*R'') , 

a{p) = l,b{p') = 1 , b^Ta^ = OL.^L2{h'^)} . 

If WF'^(T) n {W xV) mWxV, there exists a family of tempered operators T{h) : — )■ L^, 
such that T and T are microlocally equivalent on V, while T is Os^'^s^{h°°) outside V, 
that is 

for all a G S'(T*M'^) such that suppa 0^ = 0. This family, which is unique modulo 
Oy^y{h°°), is an extension of the microlocally defined T{h), see [15], Chapter 10]. 

3.3. Local /i-Fourier integral operators. We first present a a class of globally defined 
/i- Fourier integral operators following |12] and [151 Chapter 10]. This global definition will 
then be used to define Fourier integral operators microlocally. 

Let {A{t))t(z[^i^i] be a smooth family of selfadjoint pseudodifferential operators, 
Vt G [-1, 1], A{t) = Opn«(t)) , ait) G S,iiT*R^; R) , 

where the dependence on t is smooth, and WFh{A{t)) C ^2 d T*R'^, in the sense of f l3.5p . 
We then define a family of operators 

, , U(t) : L2(M'^) -> l2(M'^) , 

3.10 

^ ^ /iAf/(t) + f/(t)^(t) = 0. f/(0) = Jrf. 

An example is given by A(t) = A = a^, independent of t, in which case U{t) = exp{—itA/h). 

The family ([/(t))(g[_i^i] is an example of a family of unitary /i-Fourier integral oper- 
ators, associated to the family of canonical transformations K{t) generated by the (time- 
dependent) Hamilton vector fields Hao{t)- Here the real valued function ao(t) is the principal 
symbol of A{t) (see (13.21] ). and the canonical transformations K{t) are defined through 

|K(t)(p) = («:(t)),(i7,„(,)(p)), «:(0)(p) = p, pGT*R^ 

If f/ = f/(l), say, and the graph of k(1) is denoted by C, we conform to the usual notation 
and write 

?7g/°(M'^xM'^;C"), where C = {{x.^-y, -r,) : {x,i) = K{y,r,)} . 
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Here the twisted graph C is a Lagrangian submanifold of T*{W^ x R''). 

In words, [/ is a unitary /i-Fourier integral operator associated to the canonical graph 
C (or the symplectomorphism k(1) defined by this graph). Locally all unitary /i- Fourier 
integral operators associated to canonical graphs are of the form f/(l), since each local 
canonical transformation with a fixed point can be deformed to the identity, see [1^ Lemma 
3.2]. For any x ^ S{T*W^), the operator f/(l) x"', with x e S{T*W^) is still a (nonunitary) 
/i-Fourier integral operator associated with C . The class formed by these operators, which 
are said to "quantize" the symplectomorphism k = depends only on k, and not on the 
deformation path from the identity to k. This can be seen from the Egorov characterization 
of Fourier integral operators - see [121 Lemma 2] or [T5| §10.2]. 

Let us assume that a symplectomorphism k is defined only near the origin, which is 
a fixed point. It is always possible to locally deform /t to the identity, that is construct 
a family of symplectomorphisms K{t) on T*M'^, such that k(1) coincides with k in some 
neighbourhood V of the origin [HI Lemma 3.2]. If we apply the above construction to get 
the unitary operator f/(l), and use a cutoff x ^ S{T*W^), suppx (s V, then the operator 
U{l)x^ is an /i-Fourier integral operator associated with the local symplectomorphism 
n \ V. Furthermore, if there exists a neighbourhood V (£ V such that x\ V = I5 then 
U{l)x^ is microlocally unitary inside V. 

For an open set \/ d M"^ and K, a symplectomorphism defined in a neighbourhood V of 
V, we say that a tempered operator T satisfying 

is a micrololocally defined unitary /i-Fourier integral operator in V, if any point p & V has 
a neighbourhood Vp G V such that 

T : HiVp) H{K{Vp)) 

is equivalent to a unitary /i- Fourier integral operator associated with V^, as defined by 
the above procedure. The microlocally defined operators can also be obtained by oscillatory 
integral constructions — see for instance [211 §4.1] for a brief self-contained presentation. 

An example which will be used in §4. H is given by the standard conjugation result, see [l2l 
Proposition 3.5] or [T51 Chapter 10] for self-contained proofs. Suppose that P G \I/™'''(M'^) 
is a semi-classical real principal type operator, namely its principal symbol p = cThiP) is 
real, independent of h, and the Hamilton fiow it generates has no fixed point at energy 
zero: p = =^ dp 7^ 0. Then for any po ^ P~^(0), there exists a canonical transformation, 
K, mapping V = neigh((0, 0), r*M'^) to k{V) = neigh(po, T^M*^), with k(0,0) = po and 

poK{p)=i^{p) peV, 

and a unitary microlocal /i-Fourier integral operator U : H{V) — )■ H{n{V)) associated to 
K, such that 

U*PU = hD^^ : H{V) H{V) . 
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While C,n is the (classical) normal form for the Hamiltonian p in V, the operator hD^^ is 
the quantum normal form for P, microlocally in V. 

The definition of /i-Fourier integral operators can be generalized to graphs C associated 
with certain relations between phase spaces of possibly different dimensions. Namely, if a 
relation C C T*R'^ x T*R'' is such that its twist 

C' = {{x,^;y,-rj); {x,^;y, -rj') E C} 

is a Lagrangian submanifold of T*(M'^ x M'^), then one can associate with this relation 
(microlocally in some neighbourhood) a family of /i- Fourier integral operators T : L^(]R'^) i— 
L^iR'^) Definition 4.2]. This class of operators is denoted by IliR'^ x M*^; C), with r G M. 
The important property of these operators is that their composition is still a Fourier integral 
operator associated with the composed relations. 

3.4. Complex scaling. We briefly recall the complex scaling method of Aguilar-Combes 
[1] - see |ll],|39], and references given there. In most of this section, this scaling is inde- 
pendent of h, and allows to obtain the resonances (in a certain sector) for all operators 
P{h), h G (0, 1], where P{h) satisfies the assumptions of §2.11 

For any < ^ < ^-nd -R > 0, we define Tg^ji C C" to be a totally real deformation of 
M", with the following properties: 

TgnBc40,R) = B^40,R), 

(3.11) Te n C" \ fic"(0, 2R) = e''^M" n C" \ 5cn(0, 2R) , 

Te = {x + zfeA^) ■ a: G M"} , feA^) = Oa{e) . 

If R is large enough, the coefficients of P continue analytically outside of -8(0, R), and we 
can define a dilated operator: 

where P is the holomorphic continuation of the operator P, and u is an almost analytic 
extension of m G C^^Te^ from the totally real submanifold Te^ji to neigh(re,i?,, C"). 

The operator Pgji — z is a Fredholm operator for 26 > a.Tg{z + 1) > —26. That means 
that the resolvent, [Pg^ji — z)~^ , is meromorphic in that region, the spectrum of Pg^R in that 
region is independent of 6 and R, and consists of the quantum resonances of P. 

To simplify notations we identify ^ with M" using the map, Sg^R : Tq r — j- M", 

(3.12) Te,R3x\ — ^RexGM", 

and using this identification, consider Pg R as an operator on M", defined by {SgliyPg^RSg r 
(here S* means the pullback through S) We note that this identificaton satisfies 

C""^ ||m||l2{R") < ll'S'e,R'"||L2(re,j{) < C* ||w||l2{R") , 
with C independent of ^ if < ^ < ^o- 
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Figure 6. The complex scaling in the 2;-plane used in this paper. 

The identification of the eigenvalues of with the poles of the meromorphic continu- 
ation of 

from {Im;z > 0} to D{0, sin(26')), and in fact, the existence of such a continuation, follows 
from the following formula (implicit in [39], and discussed in |15]): if x ^ C^(M'^), supp x 
B{0,R), then 

(3.13) x{Pe,R - z)-\ = x{P - z)-\ . 

This is initially valid for Im 2 > so that the right hand side is well defined, and then by 
analytic continuation in the region where the left hand side is meromorphic. The reason 
for the Fredholm property of {Pe,R — z) in D(0, sin(2^)) comes from the properties of the 
principal symbol of Pq^r - see Fig. [61 Here for convenience, and for applications to our 
setting, we consider Pq r as an operator on L'^iW^) using the identification Sq^r above. Its 
principal symbol is given by 

(3.14) p,,^(x, = V{x + z/e,i?(x), [(1 + idfeAx))']-^^) , (x, G T*W' , 

where the complex arguments are allowed due to the analyticity of p{x,S,) outside of a 
compact set — see §2.11 We have the following properties 

Repe,ii(a;, = p{x, + 0{e^) {^f , 

Impe,H(x,0 = -4p(^,0[#e,K(x)*e] + 4p(a;,0[/e,K(a;)] + 0(^^ 
This implies, for R large enough, 

(3.16) b(x,OI<<^, \x\>2R ^ \mpe,R{^,i)<-Ce. 



For our future aims, it will prove convenient to actually let the angle 9 explicitly depend 
on h: as long as ^ > c/i log(l//i), the estimates above guarantee the Fredholm property 
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of {Pb,r — z) for z G D{0,6 /C), by providing approximate inverses near infinity. We will 
indeed take 9 of the order of h\og{l/h), see f l3.3ip . 

3.5. Microlocally deformed spaces. Microlocal deformations using exponential weights 
have played an important role in the theory of resonances since [21]. Here we take an 
intermediate point of view [26| B3] by combining compactly supported weights with the 
complex scaling described above. We should stress however that the full power of [2T] 
would allow more general behaviours of p{x,C,) at infinity, for instance potentials growing 
in some directions at infinity. 

Let us consider an independent real valued function Gq G C^(T*M'^; M), and rescale it 
in an /i-dependent way: 

(3.17) G{x, = Mh \og{l/h)Go{x, , M > fixed. 

For A G ^'™''°(]R'^), we consider the conjugated operator 

^"^ (-1)^ fl V 
= ^ ( 7^^^G»(x,hD)J A + Rl, 



(3.18) 



where 



Rl = I'e-^o-i^^m Qadc.(.,,,,)y Ae*«"(^'''^)rft. 



The semiclassical calculus of pseudodifferential operators [121 Chapter 7], [151 Chapter 4, 
Appendix D.2] and (13.171) show that 

{^^dG^i.,hD^ A = {M\og{l/h)Y{^dG^^.,hD)YA G (M/ilog(l//i))^vl>-~'0(M'^) , > . 

Since ||Gq ||i2_j.i2 < Cq, functional calculus of bounded self-adjoint operators shows that 

||exp(±tG'"'(x,/iD))|| < 

so we obtain the bound, 

with (5 > arbitrary small. Applying this bound, we may write (I3.18p as 

(3.19) ^-G-i.,hD)/hj^^G-i.,hD)/h ^ ^ Hll! QadG^^.,/.!,))'^ G Vl/-.0(M'^). 

In turn, this expansion, combined with Beals's characterization of pseudodifferential oper- 
ators (13. 4p . implies that the exponentiated weight is a pseudodifferential operator: 

(3.20) exp(G'"(x, hD)/h) G *°'^°*^(M'^) , V(5 > . 
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Using the weight function G, we can now define our weighted spaces. Let if^(M'^) be the 
semiclassical Sobolev spaces defined in (13. 3p . We put 

(3.21) H'ai^'') = e^"{--^^)Ai/,^(M'^) , He-^^'^-'^^^/^H^. , 
and 

As a vector space, if^(M^) is identical with Hj;{M.'^), but the Hilbert norms are different. 
In the case of L^, that is of /c = 0, we simply put Hq = Hq- 

The mapping properties of P = p'^{x,hD) on Hg(M.'^) are equivalent with those of 
d^t ^-G-/hp^G-/h L\R''), which are governed by the properties of the (full) symbol 
Pg of Pg- formula (I3.19P shows that 

(3.22) pG=p-iHpG + OihHog'il/h)). 

At this moment it is convenient to introduce a notion of leading symbol, which is adapted 
to the study of conjugated operators such as Pg- For a given Q G S(T*M'^), we say that 
q e S{T*R'^) is a leading symbol of Q'"{x, hD), if 

(3.23) V7 e (0, 1) , Va, /3 G N^ h-^d^d^^iQ - q) = ^./^((e)-'^') , 

that is, {Q — q) G S^~'^{T*W^) for any 7 G (0, 1). This property is obviously an equivalence 
relation inside S{T*W^), which is weaker than the equivalence relation defining the principal 
symbol map on (see §3.11) . In particular, terms of the size h\og{l/h) are "invisible" to 
the leading symbol. For example, the leading symbols of pg and p are the same. If we can 
find q independent of h, then it is unique. 

For future use we record the following: 

Lemma 3.1. Suppose 

g"'(x, hD) : Hg{W^) Hg{R'^) , Q G S{T*R'^) , 

is self-adjoint (with respect to the Hilbert norm on Hg)- Then this operator admits a real 
leading symbol. Conversely, if q & S{T*W^) is real, then there exists Q G S{T*W^) with 
leading symbol q, such that Q'^{x, HD) is self-adjoint on Hg{R'^)- 



Proof. This follows from noting that 

g;:^ = e-«"'/^Q'"(x,/iD)e^"/\ 

has the same leading symbol as Q'^{x, hD), and that self-adjointness of on Hg is 
equivalent to self-adjointness of Qg on L^: the definition of Hg in (13.211) (the case of 
A; = 0) gives 

□ 
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The weighted spaces can also be microlocahzed in the sense of §3.2t for V d T*M.'^, we 
define the space 

(3 24) HG{V) = {u = u{h)eHG{M''), : 3^ > , V/iG(0,1], < ^ 

BxeC^^iV), x'"u = u + Oy{h°^)}. 

In other words, Hg{V) = e'^'"^^''^^''^^H(V). This definition depends only on the values of 
the weight G in the open set V. 

For future reference we state the following 

Lemma 3.2. Suppose T : H{V) — )■ H{k,{V)) is an h-Fourier integral operator associated 
to a symplectomorphism k (in the sense of H3.3\) . and is asymptotically uniformly bounded 
(m the sense of i^). Take Go G {neigh{K{V))) , G = Mhlog{l/h)Go. 
Then the operator 

(3.25) T : H,^g{V) ^ Hg{^{V)) 

is also asymptotically uniformly bounded with respect to the deformed norms. 

Proof. Since the statement is microlocal we can assume that V is small enough, so that 
T = TqA in V, where Tq is unitary on L^(]R'^) and A G "^h- As in the proof of Lemma [3.11 
the boundedness of (13.251) is equivalent to considering the boundedness of 

where 

Because of (13.191) . we have uniform boundedness of Ak*g on L^. Unitarity of Tq means that 
it is sufficient to show the uniform boundedness of 

rp~l^~G'"{x,hD)/hrp^^{K''Gy"{x,hD)/h ^ ^- A'l log{l / h){Tg^ G}^ {x ,hD)To) logil / h){K* Gq)"^ (x ,hD) 

on L^. Egorov's theorem (see [151 §10-2]) shows that 

T,'G^{x, hD)To = G,{x, hD) , - k*Go e . 

Since [G^,^*^^] = h'^B B G *~°°'°(M^), the Baker- Campbell-Hausdorff formula for 
bounded operators showqj that 

rp-l^-G-^{x,hD)/hrp^^{K*G)^{x,hD)/h ^ ^- M \og{l / h)G'^ {x ,hD) ^M\og{l/h){K*GoT{x,hD) 

^ ^M\og{l/h){-G^{x,hD)+K'Gor{x,hD))+Oj^2^L^{log{l/hfh^) 

= expOL2^i2(/llog(l//l)) 
= ld + OL^-^L-ih\og{l/h)). 



^Alternatively, we can compare exp(Mlog(l//i)G'^) with (exp{Mlog{l/h)GK,))^ and use product for- 
mula for pseudodifFerential operators - see [331 Appendix] or (TSl Section 8.2]. 
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This proves uniform bounded of globally defined operators TqA, and the asymptotic uni- 
formly boundedness in the sense of fl3.8p of T on spaces of microlocally localized func- 
tions. □ 

3.6. Escape function away from the trapped set. In this section we recall the con- 
struction of the specific weight function G which, up to some further small modifications, 
will be used to prove Theorems [T] and El 

Let Ke C p~^{E) be the trapped set on the .E-energy surface, see (11.41) . and define 

(3.26) K = Ks = U Ke . 

\E\<5 

The construction of the weight function is based on the following result of Appendix]: 
for any open neighbourhoods U,V of K, U C V, there exists Gi E C°°{T*X), such that 

(3.27) Gi\u=0, HpGi>0, HpGi\p-i(^[_25,25])^ C , HpGi\p-i(^[_s,5])\v^ ■ 

These properties mean that Gi is an escape function: it increases along the flow, and strictly 
increases along the flow on p~^{[—6,6]) away from K (as specified by the neighbourhood 
V). Furthermore, HpG is bounded in a neighbourhood of p~^{0). 

Since such a function Gi is necessarily of unbounded support, we need to modify it 
to be able to use iJ^-norms defined in §3.51 (otherwise methods of [21] could be used 
and that alternative would allow more general behaviours at infinity, for instance a wide 
class of polynomial potentials). For that we follow [lH §§4.1,4.2,7.3] and [311 §6.1]: Gi is 
modified to a compactly supported G2 in a way allowing complex scaling estimates fl3.16p 
to compensate for the wrong sign of HpG2- Specifically, [31] Lemma 6.1] states that for 
any large i? > and 5q G (0, 1/2) we can construct G2 with the following properties: 
G2 e C^{T*X) and 

HpG2 > on T^(o 3R)X, 

(3.28) HpG2 > 1 on T*^o,3R)X n ip~\[-6, S]) \ V), 
HpG2 > -5o on T*X. 

Let 

G = Mh log(l /h)G2 , with M > a fixed constant. 
Then, in the notations of §3.5^ we will be interested in the complex-scaled operator 

Pg,R : (^n) _^ ^^(^n) ^ 

for a scaling angle depending on h: 

(3.29) e = e{h) = Mih\og{l/h), Ml > fixed. 
Inserting the above estimates in (13.221) . we get 

(3.30) I Repe,RM\ < ^2 , Re p ^ V , ^ Impe^nM < "^M > 
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provided that we choose [311 §6-1] 

(3.31) ^ > Ml > for some C > , 

Assuming that the constant Mq appearing in the statement of Theorem [1] satisfies 

< Mo < Ml 

for (5 > and h > small enough, the rectangle 71{S, Mq, h) is contained in the uncovered 
region in Fig. [6l hence the scaling by the angle fl3.29p gives us access to the resonance 
spectrum in the rectangle 7l{6, Mq, h). In §5.31 we will need to further adjust Mq with 
respect to Mi. 

3.7. Grushin problems. In this section we recall some linear algebra facts related to the 
Schur complement formula, which are at the origin of the Grushin method we will use to 
analyze the operator Pq^r. 

For any invertible square matrix decomposed into 4 blocks, we have 

P21 P22) - W 922 J ^ - "^'^ - ' 

provided that exists (which implies that ^22, and hence pn, are square matrices). We 
have the analogous formula for g^2^: 

(h2 = P22 - P2lPuPl2 ■ 



Pu P12 
P21 P22 



One way to see these simple facts is to apply gaussian elimination to 

V 

so that, if pii is invertible, we have an upper-lower triangular factorization: 
(3.32) V -- 



Pu 0\ A P11P12 
P21 1 / [0 P22 - P21P11P12 

The formula for the inverse of pn leads to the construction of effective Hamiltonians for 
operators (quantum Hamiltonians) P : "Hi — > We first search for auxiliary spaces 'H± 
and operators R± for which the matrix of operators 



P-z i?„ 




: Hi®H^ — >H2®H 



+ 7 



is invertible for z running in some domain of C. Such a matrix is called a Grushin problem, 
and when invertible the problem is said to be well posed. 
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When successful this procedure reduces the spectral problem for P to a nonlinear spectral 
problem of lower dimension. Indeed, if dim'H. = dim'H+ < oo, we write 



P-z R\ ^ ( Eiz) Ejz) 



and the invertibility of {P — z) : Hi — ?■ ^2 is equivalent to the invertibility of the finite 

dimensional matrix E_j^{z). The zeros of detE ^{z) coincide with the eigenvalues of P 

(even when P is not self-adjoint) because of the following formula: 

(3.33) ii hp - wY^dw = -ti J) E^+{wY^E'^^{w) dw , 

J z J z 

valid when the integral on the left hand side is of trace class - see [HJ Proposition 4.1] 
or verify it using the factorization (13.321) . Here §^ denotes an integral over a small circle 
centered at z. The above formula shows that dimker(P — z) = dimkeri? ^{z). 

The matrix E 1.(2) is often called an effective Hamiltonian for the original Hamiltonian 

P - see [H] for a review of this formalism and many examples. In the physics literature, 
this reduction is usually called the Feshbach method. 

We illustrate the use of Grushin problems with a simple lemma which will be useful later 

in 

Lemma 3.3. Suppose that 



R+ 

where Hj and 'H± are Banach spaces. If P~^ : 1-12 Hi exists then 

V is a Fredholm operator <^==^ i?+P~^i?_ : — )■ 7/+ is a Fredholm operator ^ 

and 

indP = indP+P"^P_ . 

Proof. We apply the factorization (I3.32p with pn = P, pu = R-, P21 = R+, P22 = 0. Since 
the first factor is invertible we only need to check the the Fredhold property and the index 
of the second factor: 

^0 -R+P-^R. 

and the lemma is immediate. □ 



4. A MICROLOCAL GRUSHIN PROBLEM 



In this section we recall and extend the analysis of [42] to treat a Poincare section S C 
p~^(0) for a flow satisfying the assumptions in §2.21 In [12] a Poincare section associated 
to a single closed orbit was considered. The results presented here are purely microlocal in 
the sense of §3.2[ first near a given component of the section, then near the trapped set 
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Kq. In this section P is the original differential operator, but it could be replaced by its 
complex scaled version Pe^R, since the complex deformation described in §3.41 takes place 
far away from Kq. Also, when no confusion is likely to occur, we will often denote the Weyl 
quantization of ^ symbol x ^ S{T*M.'^) by the same letter: x = X^- 

4.1. Microlocal study near S^. First we focus on a single component of the Poincare 
section, for some arbitrary k G {1, A^}. Most of the time we will then drop the subscript 
k. Our aim is to construct a microlocal Grushin problem for the operator 

near S = S^, where | Rez| < 5, | Imz| < Moh\og{l/h), and 6 will be chosen small enough 
so that the flow on ^^\kb.oz ^ small perturbation of ^^\ko- 

4.1.1. A normal form near E^. Using the assumption (12. lip and a version of Darboux's 
theorem (see for instance [231 Theorem 21.2.3]), we may extend the map = k : — )■ 
to a canonical transformation kj. defined in a neighbourhood of in T*M"', 

Qk = {{x,0 e T*W, e Sfc, \xn\ < e, l^nl < S} , 

such that 

(4.1) Kfc(x',0,^',0) = Kfc(x',^') e Sfc , poKfc = ^„. 

We call Qk = '^fc(^fc) the neighbourhood of in T*X in the range oiUk. The "width along 
the flow" e > is taken small enough, so that the sets {^2^, k = 1, . . . , A^} are mutually 
disjoint, and it takes at least a time 20e for a point to travel between any and its 
successors. 

The symplectic maps allow us to extend the Poincare section S to the neighbouring 
energy layers p~^{z), z G [—5,5]. Let us call 

Then, if 5 > is taken small enough, for z G [—5, 5\ the hypersurfaces 

Sfc(2) = /€fc,.(Sfc) = {^(x',0;e',^), (x',e') e Sfc} 
are still transversal to the flow in p~^(z). Using this extension we may continuously deform 

the departure sets Dj^ into Djk{z) =^ i^k,z{Djk) C Sfc(2;), and by consequence the tubes 
Tjk into tubes Tjk{z) C ^"^(2;) through a direct generalization of (12.121) . The tube Tjk{z) 
intersects ^j{z) on the arrival set Ajk{z) C notice that for 2 7^ 0, the latter is in 

general different from Kjz{Ajk) (equivalently Aj^^z) = Kj^l{Aji^{z)) is generally different 

from Ajk{0)). These tubes induce a Poincare map Fjk^z bijectively relating Djk{z) with 
Ajkiz). 
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The following Lemma, announced at the end of §2.3. ![ shows that for \z\ small enough 
the interesting dynamics still takes place inside these tubes: the trapped set is stable with 
respect to variations of the energy. 

Lemma 4.1. Provided 5 > Q is small enough, for any z G [—5,5] the trapped set <e 
UjkTjk{,z). 

As a consequence, in this energy range the Poincare map associated with ^{z) fully 
describes the dynamics on Kz- 

Proof. From our assumption in §2.H there exists a ball B{0,R) (the "interaction region") 
such that, for any E € [—1/2, 1/2], the trapped set Ke must be contained inside T^^q^^X. 
If R is large enough, any point p G p~^{z) \ T^^-q^-jX, z ^ 0, will "escape fast" in the past 
or in the future, because the Hamilton vector field is close to the one corresponding to 
free motion, 2 J2j ^j^xj- Hence we only need to study the behaviour of points in p~^{z) fl 

Let us define the escape time from the interaction region T^^-q j^^X: for any p G T^^q 

W(p) = inf{t > 0, max(|7r,<l>*(p)|, |7r,.<l>-*(p)|) > R} , 

For any £'g[— 1/2,1/2], the trapped set Ke can be defined as the set of points in p~^{E) 
for which tesdp) = oo. Let us consider the neighbourhood of Kq formed by the interior 
of the union of tubes, (UTjfc)°. By compactness, the escape time is bounded from above 
outside this neighbourhood, that is in p~^{0) fl T^^^ j^-^X \ (UTjfc)°, by some finite ti > 0. 
By continuity of the flow $*, for S > small enough, the escape time in the deformed 
neighbourhood p~^{z) fl T^^-q j^-^X \ {L\Tik{z))° will still be bounded from above by 2ti: this 
proves that d L\Tik{z). □ 

A direct consequence is that the reduced trapped sets Tj {z) = S(z) n Kz are contained 
inside Dj{z). 

For any set S{z) depending on the energy in the interval z G [—6, 6], we use the notation 

(4.2) S = [j S{z) . 

\z\<5 

We will extend the notation to complex values of the parameter z G 7^(5, Mq, h), identifying 
S{z) with S'(Re2;). 

4.1.2. Microlocal solutions near E. Let us now restrict ourselves to the neighbourhood of 
Sfc, and drop the index k. The canonical transformation k can be locally quantized using 
the procedure reviewed in §3.31 resulting in a microlocally defined unitary Fourier integral 
operator 

(4.3) U : H{U) — > H{n) , U* PU = HD^^ , microlocally in Vt. 
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For z G 7^(5, Mq, h\ we consider the microlocal Poisson operator 

(4.4) K(^) : ^ lL(M") , %{z)v^\{x\xn) = e*""^/%+(x') , 
which obviously satisfies the equation {hDx„ — z) K(2;) Vj^ = 0. 

For microlocally concentrated in a compact set, the wavefront set of K{z) Vj^ is not 
locahzed in the fiow direction. On the other hand, the Fourier integral operator U is well- 
defined and unitary only from Vt to VL. Therefore, we use a smooth cutoff function x^i 
= 1 in f2, = outside VL' a small open neighbourhood of Vt (say, such that < 2e 
inside fi'), and define the Poisson operator 

K{z)'^ xlUK{z) : H{f.) ^ H{n') . 

This operator maps any state G H{Y?) C L^(M"~^), to a microlocal solution of the 
equation (P — z)u = in f2, with u G H{Q'). As we will see below, the converse holds: 
each microlocal solution in Q is parametrized by a function G i^(S). 

In a sense, the solution u = K{z)v^ is an extension along the fiow of the transverse data 
v+. More precisely, K{z) is a microlocally defined Fourier integral operator associated with 
the graph 

(4.5) C_ = {(^(a;',x„,e',Rez);a;',0, < e} C T*(X x M^-^) . 

Equivalently, this relation associates to each point {x',^') G S a short trajectory segment 
through the point /t(x', 0; C,', Rez) G S(Re2;). We use the notation C_ since this relation is 
associated with the operator R_ defined in fl4.13p below. 

Back to the normal form hD^^, let us consider a smoothed out step function, 

Xo e C~(M^J, XoM = for x„ < -e/2, xoM = 1 for x„ > e/2 . 

We notice that the commutator {i/h)[hDx^,xo] = x'oi^n) is localized in the region of the 
step and integrates to 1: this implies the normalization property 

(4.6) {{i/h)[hD^„,xo]^{z)v+,K{z)v+) = \\v+\\l2^^„-i) , 

where (•, •) is the usual Hermitian inner product on L^(M'"). Notice that the right hand 
side is independent of the precise choice of Xo- 

We now bring this expression to the neighbourhood of S through the Fourier integral 
operator XqU. This implies that the Poisson operator K{z) satisfies: 

(4.7) {{z/h)[P,xnK{z)v+,K{z)v+) = \\v+f for any v+ G iJ(S) . 

Here the symbol x is such that x^ = U Xo U* inside Q, so x is equal to before $~'^(S) 
and equal to 1 after $^(S) (in the following we will often use this time- like terminology 
referring to the fiow $*). In (14. 7p . we are only concerned with [P, x^] microlocally near Q, 
since the operator Xn^ i^ microlocalized in Q' x Q'. Hence, at this stage we can ignore the 
properties of the symbol x outside Q'. 
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The expression (14. 7p can be written 

(4.8) K{zy [{t/h)P, x^Kiz) = Id : i/(S) ^ if(S) . 

Fixing a function x with properties described after (14. 7p and writing x = Xf (where / is 
for forward), we define the operator 

(4.9) R^{z) K{-zy [{z/h)P, Xf] = K{zy U*Xn [{^/h)P, Xf] 

(from here on we denote x = X^ similar expressions). This operator "projects" any 
u G H{Q) to a certain transversal function v+ G -f^(S). But it is important to notice that 
R+{z) is also well-defined on states u microlocalized in a small neighbourhood of the full 
trapped set K: the operator Xn [{i/h)P, Xf\ cuts off the components of u outside VL. Hence, 
we may write 

R+{z) ■ H{neigh{K)) H{T,) . 

The equation (14.81) shows that this projection is compatible with the above extension of 
the transversal function: 

(4.10) R+{z)K{z) = Id : i/(S) ^ if(S) . 

This shows that transversal functions v + G -ff (S) and microfocal solutions to (P — z)u = 
are bijectively related. Since |Im^| < Moh\og{l/h) and \xn\ < 2e inside fl' (resp. \xn\ < e 
inside fl), we have the bounds 

Just as K{zy, R+{z) is a microlocally defined Fourier integral operator associated with the 
relation 

(4.11) C+ = {x,^';{K{x\xn,^',Rezy, (x', x„, e', Re 2) G fi} C T*{W'-'^ x X), 

namely the inverse of C_ given in (14.51) . In words, this relation consists of taking any 
p G np~^{Rez) and projecting it along the flow on the section T,{z). 

We now select a second cutoff function Xb with properties similar with Xf, and satisfying 
also the nesting property 

(4.12) Xb = 1 ill a neighbourhood of suppx/ • 
With this new cutoff, we define the operator 

(4.13) R4z)u^ = [{i/h)P, Xb] K{z) : Hijl) ^ H{Sl) . 

Starting from a transversal data m_ G H{Tj), this operator creates a microlocal solution in 
Vt and truncates by applying a pseudodifferential operator with symbol HpXb- Like K{z), it 
is a microlocally defined Fourier integral operator associated with the graph C_. its norm 
is bounded by \\R^{z)\\l2^l^ = 0{h-'^^^). 
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4.1.3. Solving a Grushin problem. We are now equipped to define our microlocal Grushin 
problem in Q. Given v G H{Q), v+ G -f^(S), we want to solve the system 



(4.14) 



{i/h){P - z)u + R_{z)u_ =v, 



with u G L'^{X) a forward solution, and u_ G -f^(S). 

Let us show how to solve this problem. First let u be the forward solution of {i/h){P — 
z)u = V, microlocally in Q. That solution can be obtained using the Fourier integral oper- 
ator U in (14. 3 p and the easy solution for hD^^. We can also proceed using the propagator 
to define a forward parametrix: 

(4.15) u = E{z)v, E{z)''^ f e-^'^''-'^'^ dt. 

Jo 

The time T is such that $^(f2) fl = (from the above assumption on the separation 
between the Qk we may take T = 5e). By using the model operator hDx„, one checks that 
the parametrix E{z) transports the wavefront set of v as follows: 

(4.16) WFh{E{z)v) cWFh{v)U^'^{WFh{v))U [j ¥ {WF h{v) f] p-\Re z)) . 

0<t<T 

In general, u does not satisfy R+{z)u = so we need to correct it. For this aim, we solve 
the system 



(4.17) 

through the Ansatz 
(4.18) 



(i/h) (P - z)u + {z)u^ = , 
R+{z)u = — R^{z)u 



= -f+ + R+{z)u , 
= -XhK{z) M_ . 

Indeed, the property (P — z) K{z) = ensures that {i/h){P — z)u = —R_{z)u_. We then 
obtain the identities 

R+{z)u = -K[-zy [{i/h)P, Xf] Xh K{z) u. 
^-K{-zy[{z/h)P,Xf]K{z)u^ 

The second identity uses the nesting assumption {HpXf)Xb = HpXf, and the last one 
results from (14.81) . This shows that the Ansatz (I4.18P solves the system (14.170 . Finally, 
{u = u+u, U-) solves (I4.14P microlocally in f2xS, for v G H{Q) and G H{T,) respectively. 
Furthermore, these solutions satisfy the norm estimate 

(4.19) ll^ll + < /i"5^^«^(||t;|| + 
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The form of the microlocal construction in this section is an important preparation for the 
construction of our Grushin problem in the next section. In itself, it only states that, for v 
microlocalized near S, {i/h){P — z)u = v can be solved microlocally near S in the forward 
direction. 



4.2. Microlocal solution near K. We will now extend the construction of the Grushin 
problem near each S^, described in §4.11 to obtain a microlocal Grushin problem near the 
full trapped set K. This will be achieved by relating the construction near to the one 
near the successor sections Sj. We now need to restore all indices k G {1, . . . , A^} in our 
notations. 



4.2.1. Setting up the Grushin problem. We recall that H(T,k) C L'^{W^~^) is the space 
of functions microlocally concentrated in (see (13. 6p ). For u G L'^{X) microlocally 
concentrated in neigh(i^', T*X), we define 



(4.20) R+{z)u = {Rl{z)u, R+{z)u) G i/(Si) x ... x i7(S 



N 



where each R\{z) : if(neigh(i^)) — )■ if(Efc) was defined in §4.11 using a cutoff X/ ^ 
C'^{T*X) realizing a smoothed-out step from to 1 along the flow near S^. 

Similarly, we define 



R.{z) : F(Ei) X ... X i7(Ejv) ^ H {u'^^,^!^) , 

(4.21) J!, . 

R_{z)u^ = ^Rl{z)ul, u_ = (ui,...,u^^ 



Each R^i{z) was defined in (14.131) in terms of a cutoff function x\ ^ C'^{T*X) which also 
changes from to 1 along the flow near S^, and does so before X/- Below we will impose 
more restrictions on the cutoffs xt- 

With these choices, we now consider the microlocal Grushin problem 



(4.22) 



{i/h){P - z)u + R^{z)u^ =v, 
R+{z)u = f+ 



The aim of this section is to construct a solution {u, uJ) microlocally concentrated in a 
small neighbourhood of 

Ko X K-^^iTi) X ... X «:-^(7]v) , 
provided (f,f+) is concentrated in a sufficiently small neighbourhood of the same set. 

To achieve this aim we need to put more constraints on the cutoffs xt- assume that 
each Xb ^ C'^{T*X) is supported near the direct outflow of 71-. To give a precise condition. 
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let us slightly modify the energy-thick tubes Tj^ (see fl2.12p . (14. 2p ) by removing or adding 
some parts near their ends: 



jrs.s, dcf ^^t^p^ . p ^ -^^^^ 2e<t< t+{p) + s,2e}, = ± . 

With this definition, the short tubes TJJT do not intersect the neighbourhoods fi^, fij, while 
the long tubes T^^'*' intersect both (see Fig. [7]). 
We then assume that 

(4.23) Xtip) = l for pG (J f-', 

and suppXb is contained in a small neighbourhood of that set. Furthermore, we want 
the cutoffs {Xb)k=i,...,N to form a microlocal partition of unity near Kq: there exists a 
neighbourhood Vq of K containing all long tubes: 



(4.24) VoD[jf+-^ 
and such that 

N 

(4.25) ^;^^(p)^l for peVo. 

k=l 

These conditions on Xb can be fulfilled thanks to the assumption (12.101) on the section 
S. A schematic representation of these sets and cutoffs is shown in Fig. [3 

4.2.2. Solving the homogeneous Grushin problem. Let us first solve fl4.22p when v = 0. The 
wavefront set WF/i(u^) C is mapped through k^^z to a subset of T,k{z). The microlocal 
solution Kk{z)v^, initially concentrated inside the neighbourhood can be extended 
along the flow to a larger set fi^, which intersects the successors T,j{z) of T^k{,z) and contains 
the union of tubes Ujgj+(fc) ^jfc"'' ('^^ remind that j ^ k according to assumption f l2.10p ). 
This can be done by extending the symplectomorphism k^, the associated unitary Fourier 
integral operator f/^, and replace the cutoff function by a function Xn+ supported in 

k 

the set we can then define the extended Poisson operator as: 

Kt{z) = xl^^ U, Kiz) : i7(S) ^ Hint) . 

Assuming Kfc,z(WF/i(w^)) is contained in the departure set Dk{z) C Sfc(z), the extended 
microlocal solution Kt{z)v^ is concentrated in the union of tubes Uj^j^(^k)T^f^ {z). In that 
case, we take as our Ansatz 

(4.26) u,'^xlKt{z)vX. 
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Figure 7. Schematic representation of (part of) the neighbourhoods Vi C 
Vq of Kq (resp. green shade and green dashed contour), some sections 
(thick black) and arrival sets Akj C (red). We also show the tubes Tj|^ 
connecting S2 with A12 (the dashed lines indicate the boundaries of Tj^~), 
the supports of the cutoffs Xb xj (dot-dashed line), and two trajectories 
in Kq (full lines inside Vi). 

Set 

(4.27) tmax = max{t+(p), p e UkDk{z), \Rez\ < 6} 

the maximal return time for our Poincare map. Then the above Ansatz satisfies the estimate 

(4.28) \\uk\\L^<h-'''^'--+'^\\v+\\H,j^.. 



Due to the assumption (14. 23 p . the cutoff xt effectively truncates the solution only near 
the sections ^^(2;) and T,j{z), j G J+{k), but not on the "sides" of suppxt- Hence, the 
expression 

(4.29) {z/h)iP - z)u, = [{i/h)P, xl] Kt{z) vX 
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can be decomposed into one component R'^_{z)v']^ supported near Dk{z), and other com- 
ponents supported near the arrival sets Ajk{z) C flj, due to the "step down" of xt ^^ar 
Ajk{z). The assumption (14.25P ensures that 

(4.30) [{i/h)P, xl] = -[{^lh)P, x^,] microlocally near A^k{z) , 
so the expression in fl4.29p reads 

(4.31) {tlh){P-z)u,^Bt{z)vX- J2 [m)P,xi\K^{z)vi. 

Now, for each j G J+{k) we notice that K^{z) is a solution of {P — z)u = near Ajk{z), 
so this solution can also be parametrized by some transversal data "living" on the section 
^j{z) (see the discussion before fl4.5p ). This data obviously depends linearly on f^, which 
defines the monodromy operator M.jk{z): 

(4.32) K^{z)v^^ = Kj{z) Mjk{z)v'l , microlocally near Ajk{z) . 

The operators Aijk{z) are microlocally defined from C to Ajk{z) C Sj, they are 
zero on H{D^k) for I ^ j. The identity fl4.8p provides an explicit formula: 

(4.33) M,u{z) = K,{zy [{z/h)P, x^j]K^{z) = Ri{z)K^{z) . 

Before further describing these operators, let us complete the solution of our Grushin 
problem. Combining (14.311) with (14.321) . we obtain 

(4.34) (^^/h){P-z)uk = R'l{z)v1- J2 R-i^WjkizK . 

This shows that the problem (I4.22p in the case v = and a single v'^, WF/,(t>^) C D^. is 
solved by 

u = xt Kt{z) vX, ut = -vl, uL = Mjk{z)vX, J e J+{k) . 
We now consider the Grushin problem with v = 0, = {v]_, . . . ,v^) with each mi- 
crolocalized in Dk- By linearity, this problem is solved by 

u^J2^tKUz)v1, 

(4.35) 



keJ-{j) 



From the above discussion, u is microlocalized in the neighbourhood Vq of K, while uL is 



microlocalized in Dj U Aj(z) 



Let us now come back to the monodromy operators. The expression fl4.33p shows that 
A4jk{z) is a microlocal Fourier integral operator. Since we have extended the solution 
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Kk{z) beyond Qk, the relation associated with the restriction of K^{z) on H{Djk) is a 
modification of (14. 5p . of the form 

Ci^ = {($*(k,.^,(p)); p), p g D,k, -e<t< t^ax + e} , 

such that the trajectories cross Sj. On the other hand, the relation C+ associated with 
R''^{z) is identical with (14. lip . By the composition rules, the relation associated with 
Mjk{z) is 

C^^ = {(p', p), p G Djk, P = i^ll o Fjk,^ o KkAp) = Fjk,z{p)} ■ 

This is exactly the graph of the Poincare map -Fjfc,z : Djk{z) — )■ Ajk{z), seen through the 
coordinates charts k^^z, i^j,z- 

When z is real, the identity (14.81) implies that M.jk{z) : H{Djk) — )■ H{Ajk{z)) is microlo- 
cally unitary. Also, the definition (I4.33P shows that this operator depends holomorphically 
of z in the rectangle TZ{S, Mq, h). To lowest order, the z-dependence takes the form 

M.kiz) = M.kiO) 0pnexp(z2t+//i)) (1 + 0(/ilog(l//i))) 

where G C^{W^~^] M+) is an extension of the return time associated with the map Fj^^z 
on Djk- For z G 7^(5, Mq, h), this operator satisfies the asymptotic bound 

(4-36) ll-^.;^(^)L(5,)^^(A~ (.)) = O{h-'^ot^-) 

4.2.3. Solving the inhomogeneous Grushin problem. It remains to discuss the inhomoge- 
neous problem 

(4.37) {i/h){P- z)u + R^u^ = v, 

for V microlocalized in a neighbourhood Vi of K, which satistfies 

(4.38) V^c[jfJi^. 

(each tube T~f^^ intersects fifc only near D^, see figure [7]). 

Let us first assume that v is microlocally concentrated inside a short tube Tjj^^ . We use 
the forward parametrix E{z) of {i/h){P — z) given in (I4.15P with the time 

(4.39) T = t^ax + 5e , 
and consider the Ansatz 

(4.40) u^='xtE{z)v. 

According to the transport property (I4.16p . E{z)v is microlocalized in the outflow of T~ir , 
so the cutoff Xb effectively truncates E{z)v only near Ajk{z) C flj. The partition of unity 
(I4.25P then implies that 

iz/h){P -z)u = v + [{i/h)P, xl] E{z) V = v- [{i/h)P, x^J E{z) v . 



FROM OPEN QUANTUM SYSTEMS TO OPEN QUANTUM MAPS 37 

Also, the choice of the time T ensures that E{z)v is a microlocal solution of (P — z)u = 
near Ajk{z), so 

E{z)v = Kj{z)R-'_^{z)E{z)v microlocally near Ajk{z) . 
Thus, we can solve fl4.37p by taking 

ui = lV_^{z)E{z)v, ni = 0, l^]. 

The propagation of wavefront sets given in (14.161) shows that u'_ G H{Ajk{z)), and that 
WF/i(n) C TyfT does not intersect the "step up" region of the forward cutoffs Xj, so that 
i?^(z)n = Ofor alH = l,...,iV. 

If V is microlocally concentrated in Vi fl U|t|<e$*(-Dfc), we can replace the cutoff xt 
(ICTll by 

and apply the same construction. The only notable difference is the fact that B!^{z)u may 
be a nontrivial state concentrated in [J\t\<eDj^. 

In both cases, we see that + < /i~^"^°'^*™="'"''^'^)||t>||. By linearity, the above 

procedure allows to solve (I4.37P for any v microlocalized inside the neighbourhood Vq. 

This solution produces a term which can be solved away using the procedure of 

g221 Notice that \\R+u\\ < 

We summarize the construction of our microlocal Grushin problem in the following 

Proposition 4.2. For 5 > small enough, there exist neighbourhoods of K = Ks in T*X, 
V+ and and neighbourhoods of'k'J^iTj) in Sj, V^, and Vl, j = 1, ■ ■ ■ N , such that for 
any 

{v, v+) G H{V+) X H{Vl) X ■ ■ ■ H{V^) , 

we can find 

G i/(Vl) X H{V}) X ■ ■ -//(v:^) , 

satisfying 

— (P — z)u + R-{z)u- = V , R+{z)u = v+ microlocally inV^xV^x--- . 

Here R±{z) are given by (I4.20p and (I4.2ip . Furthermore, the solutions satisfy the norm 
estimates 

\M + \\uA\<h-''''^'''--^^'\\\v\\ + \\v+\\), 
where tmax is the maximal return time defined in (14.270 . 
One possible choice for the above sets is 

\Kcz\<5 
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Proof. Take v G HiVi), and call (m, m_) the solution for the inhomogeneous problem fl4.37p . 
Then the propagation estimate fl4.16p implies that u is concentrated inside the larger neigh- 
bourhood Vo C Uj,fef++ (see (jSMD), while ui G H{Aj{z)). 

We have B!^{z)u G H{Dk) so, provided the data satisfies G Dk, the computations 
of §4.2.21 show how to solve the homogeneous problem with data — R+{z)u). That 
solves the full problem. The expressions (14.351) show that the solutions to the homogeneous 
problem (u, u^) are microlocalized, respectively, in Vq and in Dk UAk{z). □ 

Remark. The proof of the proposition shows that the neighbourhoods and are 
different. For given data {v, v'^), the solutions {u, uJ) will not in general be concentrated in 
the same small set as the initial data. This, of course, reflects the fact that a neighbourhood 
V of Kq is not invariant under the forward flow, but escapes along the unstable direction. 
In order to transform the microlocal Grushin problem described in this proposition into a 
well-posed problem, we need to take care of this escape phenomenon. This will be done 
using escape functions in order to deform the norms on the spaces L'^{X) (as described in 
§3.5p . but also on the auxiliary spaces L^(M"'~^). 

5. A WELL POSED Grushin problem 

The difficulty described in the remark at the end of ^will be resolved by modifying the 
norms on the space L'^{X) x L^(M'^~^)^, through the use of exponential weight functions 
as described in §3.51 These weight functions will be based on the construction described in 

mM 

In most of this section we will consider the scaled operator Pg^n globally, so we cannot 
replace it by P any longer. To alleviate notation, we will write this operator 

(5.1) P = Pe,R, e = M^h\og{l/h) , i?>Co, 

where Co is the constant appearing in (12. 2p . and Mi > is a constant (it will be required 
to satisfy (13.311) once we fix the weight G, and is larger than Mq appearing in Theorem [1]). 

We will first discuss the local construction near each S/j and then, as in the previous 
section, adapt it to construct a global Grushin problem. 

Our first task is still microlocal: we explain how a deformation of the norm on L'^{X) by 
a suitable weight function G can be used to deform the norms on the N auxiliary spaces 
L^(]R"~-'^), microlocally near E^. 

5.1. Exponential weights near S^. As in §4.H in this subsection we work microlocally 
in the neighbourhood Qk of one component {Q^ is the neighbourhood described in §4.ip : 
we drop the index k in our notations. Notice that the complex scaling has no effect in this 
region, so P = Pq^r. We will impose a constraint on the weight function G near E, and 
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construct a weight functions (7 on S. The construction of the local solution performed in 
§4.11 will then be studied in these deformed spaces. 

Take a function g° G C^{W-^), and use it to define Go G C°°(T*M'"), so that 

Go{x',Xr„C,U=g\x',e) inn'. 

Then, using the Fourier integral operator U given in (14. 3p . one can construct a weight 
function Go G S{T*X) such that 

= U (Go)"' U* microlocally near fi. 
Notice that Gq now depends on h through an asymptotic expansion 

(5.2) Go{h) ~ ^o,i , Goj G C^{T*X) independent of h . 

j>o 

This weight satisfies Go,o = Gq o in Q, and the invariance property 

(5.3) [P{h), G^{x, hD)] = microlocally in n . 
As in §3.51 we rescale these weight functions by 

(5.4) G = M/ilog(l//i)Go, g = Mh\og{l/h)g''. 
Still using the model hDrc^, one can easily check the intertwining property 

G"'(x, hD^; h) K{z) = K{z) g'"{x\ hD^r, h) : HiT.) H{Q') , 

Using the weights G and g we define the microlocal Hilbert spaces Hg{^') and HgiY?) by the 
method of §3.5[ We need to check that the construction of a microlocal solution performed 
in §4. 1.21 remains under control with respect to these new norms. 

Lemma 5.1. The operators 

K{z) : Hg{E) Hcin') , z G 7^(5, Mq, h) 

satisfy the analogue of (14.71) . Namely, taking a cutoff x jumping from to 1 near E as in 
then any f+ G Hg{T,) will satisfy 

(5.6) {[it/h)P, K{z) v^, K{z) v^)Ha ^Ml^. 

Proof. From the cutoff x W6 define the deformed symbol xg through 

xS(x, hD) e-^'^^^''^^)/'^ x"(a:, hD) eG™(-''^^)/^ . 
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The symbol calculus of §3.5l shows that xg also jumps from to 1 near S, so that (returning 
to the convention of using x for X^) 

([(z//.)P, xmz)v+. K{z)v+)Ha = (e^^/^(z//^)P, x]K{z)v^. e-«/^ir(z)i;+)^2 

^ {K{zy[{i/h)PG, Xg] K{z) e~^l^ e~^l^ v+)l2 
^ {K{zr[{i/h)P, Xg] K{z) e-^l^ v+, e~^/^ v+)l2 

In the second line we used (15.51) . the third line results from P = Pq, due to (15.31) . and the 
last one from (14. 7p applied to xg- D 

Equation (15. 5p shows that, for z G 7l{6, Mq, h), the operators K{z) and R±{z) defined 
respectively in (14. 9 p and (14.131) . satisfy the same norm estimates with respect to the new 
norms as for the norms: 

(5.7) II^WIlH,(£HH,(n) = (^(^-''"), 

(5.8) \\R4^)\\Hain)^H,i^) = o{h-^n , ii^-(^)ii^,(£)^^,,(n) = o{h-''n ■ 

The arguments presented in §4.11 carry over to the weighted spaces, and the microlocal 
solution to the problem (I4.14p constructed in §4.1.31 satisfies the norm estimates 

(5.9) \\u\Ua + \\uA\h, < h-''''>^{\\v\\Ha + \Mh,) . 

Given a function 6*0.0(3^,0 satisfying HpGofi = in f2, one can iteratively construct a full 
symbol Go of the form (15.20 . such that (15.30 holds. Now, the lower order terms in Go may 
change the norms only by factors (l + 0{Mhlog{l/h))^, so the same norm estimates hold 
if we replace Go by its principal symbol Go,o in the definition of the new norms. As a result, 
we get the following 

Proposition 5.2. Take E C^{T*W~^), Go{x',Xn,^',^n) = g^{x',^'), Go G C^{X) 
satisfying Gq = Gq ok in Q, and 

G = Mh\og{l/h)Go, g = Mh\og{l/h)Go. 

Then, the estimates 1 5. l- iS. 9\) hold in the spaces Hg{^), Hg{T,). 



5.2. Globally defined operators and finite rank weighted spaces. In this section 
we transform our microlocal Grushin problem into a globally defined one. This will require 
transforming all the microlocally defined operators {R±{z), A4jk{z)) into globally defined 
operators acting on L'^{X) or L^(M"'~^). Because our analysis took place near the trapped 
set Kq, we will need to restrict our auxiliary operators to some subspaces of L^(M"~^) ob- 
tained as images of some finite rank projectors. These subspaces are composed of functions 
micro localized near Kq. To show that the resulting Grushin problem is well-posed (invert- 
ible), the above construction must be performed using appropriately deformed norms on 
the spaces L^(X), L^(]R"~^), obtained by using globally defined weight functions G, gj. 
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Our first task is thus to complete the constructions of these global weights, building on 



5.2.1. Global weight functions. We will now construct global weight functions G G C^{X), 
Qj G C^{T*W^~^) (one for each section Sj). For this, we will use the construction of an 
escape function away from Kq presented in §3.6[ and modify it near the Poincare section so 
that it takes the form required in Proposition 15.21 and allows us to define auxiliary escape 
functions gj. These weight functions will allow us to to define finite ran/j realizations of the 
microlocally defined operators R±{z) and Ai{z). 

Our escape function Gq G S(T*X) is obtained through a slight modification of the weight 
G2{x,^) described in fl3.28p . The modification only takes place near the trapped set K, 
and in particular near the sections Zlj. The following lemma is easy to verify. 



Lemma 5.3. Let }j=i,...,A- be the neighbourhoods ofT^j described in i. i , and Vt"- 



be small neighbourhoods ofQj, Qj d Q'j d Qj, and let V be a small neighbourhood of Kg 
(see (13:261) ). Then there exists Gq G C^{T*X) such that 



N 

dcf 



H,Go>l on T^^,^,j,)Xnp~\[-6,6])\W, W = Vu[jn% 

i=i 

(5-10) HpGo = on Q'j, 

HpGo > on T^(o 3^)X, 
HpGo>-6o onT*X. 

Besides, using the coordinate charts : Q'j — Q'j (see ^UTTjj, we can construct Gq such 
that Gq o Q'j is independent of the energy variable C,n £ S]- 

The last assumption (local independence on is not strictly necessary, but it sim- 
plifies our construction below, making the auxiliary functions gj independent of 2; — see 
Proposition 15. 2[ 

For the set V we assume that V Vi, where Vi is the set defined in fl4.38p . As a 
consequence, there exists a set V{, with V (£ V( (£ Vi with the following property. Consider 
the the parametrix E{z) 04.151) with the time T = tjnax + 5e. Then there exists ti > such 
that, for any p G p~^{[—6,6]) \ V(, the trajectory segment {$*(p), < t < T} spends a 
time t > ti outside ofW. The main consequence of this property is a strict increase of the 
weight along the flow outside V(: 

(5.11) Vp G T*(o,2^)X n p-\[-6, 6]) \ V; , Go($^(p)) - Go(p) > h . 

(Here we use the fact that T is small enough, so that a particle of energy z ^ starting 
inside T^(^q2R) at t = will remain inside 7"^(0 3i?) up to t = T.) The set V will be further 
characterized in the next subsection. 
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From now on, we will take for weight function G = Mh log h Gq with such a function Gq, 
and use it to define a global Hilbert norm || • in fl3.2ip . As in Proposition 15. 2[ 

we define, for each j = 1, . . . , N , the auxiliary weight 

(5.12) g,{x', = Mh log(l//i) Go o (x', 0, e', 0) , (x', O e , 

and extend it to an element of C^(T*]R*^"^^)), so that the deformed Hilbert norm 

is globally well-defined. Proposition 15.21 shows that our microfocal construction near 
satisfies nice norm estimates with respect to the spaces Hg{X), Hg.. 

To see the advantages of having weights which are escape functions we state the fol- 
lowing lemma which results from applying Lemma 13.21 to the Fourier integral operator 
exp{-itP/h): 

Lemma 5.4. Suppose that pi = $*(po) for some t > 0, and that 

A = Go(pi)-Go(po) >0. 

Suppose also that Xj ^ C^{T*X), j = 0, 1, have their supports in small neighbourhoods of 
Pj 's. Then for h small enough we have 

(5.13) lle-^*^/'^ X^Ho^Ha < h''^^' , llxr e-''''/'\\Ha^Ha < h""^/' . 

5.2.2. Finite dimensional projections. We want to construct a finite dimensional subspace 
of the Hilbert space Hg.(W^~^), such that the microlocal spaces Hg.(y].) are both approxi- 
mated by it modulo 0{h°°). 

For each j = 1, . . . , A^, let S'j, Sj be two families of open sets with smooth boundaries in 
T*W-\ satisfying 

(5.14) <^ S'^ <m S, C , j = l,...,N. 

In particular, each Sj, S^ splits into disjoint components S^j d Skj C D^j. 

Once these sets are chosen, we need to adjust the set V in Lemma [573| making it thinner 
if necessary: 

Lemma 5.5. For S > small enough, there exists V = neigh(ii'5, Vi) and to > such that 
the following property holds. 

For any indices j = 1, . . . ,N, k G J+{j), any z G [—6, 6] and any point p G Dkj fl Sj 
such that its successor F^j^zip) does not belong to S'j^, then the trajectory between Kj,z{p) 
and Fkj z{K,j^z{p)) spends a time t >tQ outside ofW = VU IJj=i ^j- 

The time to is necessarily smaller than the maximal return time tmax of fl4.27p : on the 
other hand, to increases if we decrease the width ~ e of the sets Qj. See figure M for a 
sketch. Now, let 
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Figure 8. Schematic representation (inside some energy layer p~^{z)) of 
the neighbourhood V and the sets 5*^, 5*^. The departure/arrival sets D^j, A^j 
are similar to the ones appearing in figure HI The sets Sk, Sj are represented 
through their images in S^, Sj through k^. 2, k^ ^. We showe 3 trajectories 
staying inside V all the time, and one ending outside 5*^. 



Qj = Qj{x',^';h)eS{T*R'^~'), 

with leading symbol qj independent of h (the leading symbol is meant in the sense of 
f l3.23p ). We choose that leading symbol to be real and have the following properties: 

qj{p)<0, peSj, 

^^■^^^ qj{p) > , p G T*M"-^ \ , liminf gj(p) > 0. 

Lemma [3.11 shows that one can choose Qj so that 

QJ{x',hD^,) : Hg^{W~^) — >Hg^{W-^) is self-adjoint. 

Under the assumptions fIS.lSp . we know that Qj has discrete spectrum in a fixed neigh- 
bourhood of M_ when h > is small enough. Let 

(5.16) Hj = U,{Hg^{«--')) , where = {QJ{x\ hD^,)) , 

that is, II j is the spectral projection corresponding to the negative spectrum of . In 
particular, 

(5.17) \mW-^H,^ = 1, dim(H,) ~ c, h'-"^ , c, > . 
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We group together these projectors in a diagonal matrix Uh '= diag(ni, . . . , IlAr) projecting 
if3,(M"-i) X ■ • ■ iJg^(M"-i) onto H = Hi X ■ ■ ■ Hn. 

The space Hj will be equipped with the norm || • \\Hg.- For future reference we record 
the following lemma based on functional calculus of pseudodifferential operators (see for 
instance [121 Chapter 7]): 

Lemma 5.6. For any uniformly bounded family of states u = {u{h) G L'^{W^^^))h^Q, 
WYn{u) ^S, ^\\u- U,u\\h,^ = Oih^)\\u\\H,^ . 

In §5.11 we used the microlocally defined operators 
Renaming them R-^_^j^{z) (where m stands for microlocal) we now define 

(5.18) i?i(^) = : Hg{x) ^ n, . 

The estimate (15.81) together with the above Lemma shows that 

(5.19) \\R^^{z)\\Ha{x)^H, = 0{h-^'n . zen{6,Mo,h). 

The operators R\{z) are globally well-defined once we choose a specific realization of 
R'^^lz), which gives a unique definition mod 0{h°°). We have thus obtained a family 
of operators 

i?+(z) = (i?^,...,i?^) : HG{X)-^'HiX---nN. 

In turn, the operators R''_{z) are obtained by selecting a realization of the microlocally 
defined operator R-'_,^[z) on Hg.{T,j), and restricting that realization to Hj: 

(5.20) R^_{z) = Rl^^{z) li, : Hg{X) . 

Again, these operators are well defined mod 0{h°°). Putting together (15. 8p with (15.171) 
ensures that 

\\RUz)\\n,^na = 0{h~^'n ■ 

We group these operators into 

R-{z) : Hi X ■■■Un ^ Hg{X) 
(5 211 ^ 

^' ' R.{z)u^ = ^Rl{z)ul, u. = {u\,...,ui). 

i=i 
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5.3. A well posed Grushin problem. With these definitions we consider the following 
Grushin problem: 

S C'/""^';',^)' - ^' ^(^') . . e ns, Mo. k) . 

Since Pe,R{h) — z (which we will denote hj P — z for short) is a Fredholm operator, so is 
V{z), as we have only added finite dimensional spaces. For Imz > the operator (P — z) 
is invertible, so Lemma 13.31 shows that the index of V{z) is 0. Hence, in order to prove 
that V{z) is bijective it suffices to to construct an approximate right inverse and then use a 
Neumann series. The rest of this section will be devoted to the proof of this (approximate) 
right invertibility of V{z). 

5.3.1. A well-posed homogeneous problem. As before we first consider the homogeneous 
problem 

(5.23) {i/h){P - z)u + R_{z)u^ = 0, R+{z)u = v+ , 

where only one component is nonzero (we may assume that = !)• For that we 

adapt the methods of §4.2.21 We construct an approximate solution using the extended 
Poisson operator K^{z) (that operator acts on the microlocal space Hg^^ijlk), so its action 
on "Hfc is well-defined modulo (9(/;.°°)), and take 

u = xlKt{z) vX, 

where x\ is the backwards cutoff function with properties given in f l4.12p .( H.23p and f l4.25p . 
The fact that G increases along the trajectories implies that u satisfies the same norm 
bound as with the "old norms" (see f l4.28p ): 

The microlocally defined operator satisfies 

RX^^{z)u^vX + OH,Shn. R\,„Az)u = OH,^ih^), Jy^k. 

As a result, projecting the left hand side onto 'H'' has a negligible effect: 

Rliz) u = nj,{vX + 0{h^)) = vX + On.ih^) • 
Following f l4.29p we write 

(5.24) {t/h)iP - z)u ^ [{i/h)P, xl] Kt{z)v^ e Hg{X) . 

As noticed in §4.2.2[ the transport properties of K'^{z) show that u is microlocalized inside 
the union of tubes Uj^j_^_(^k)T^k^{z), so the right hand side in (15.240 splits into a component 
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concentrated near D^, and other components concentrated near the arrival sets Ajk[z), 
j G J+{k). We rewrite f l4.34p for the present data: 

(5.25) (^^/h)iP-z)u^Rt,^{z)v1- Ri^,mi^)M,,iz)v'^ . 

j&j+ik) 

Each state M.jk{z)v^ is microlocalized inside the arrival set Ajk{z) C Sj, which is not 
contained in Sj in general - see the remark at the end of §4] and Fig. |8l 

Consequently one could fear that replacing the operators R-'_ „^{z) by the truncated oper- 
ators R''_{z) would drastically modify the above right hand side. The microlocally weighted 
spaces Hqi Hg. have been constructed precisely to avoid this problem. The mechanism is a 
direct consequence of the relative properties of the sets Sj and V explained in Lemma 15.51 
Namely, a point pk G Sjk is either "good", if its image pj = Fjk^z{pk) £ S'j, or "bad", in 
which case 

(5.26) Goip,)-GoiPk)>to. 
Let us choose a cutoff 

(5.27) Xj e C^iSj), Xj = 1 on , = outside neigh(S';, S.j) . 

Since the Fourier integral operator : H{Dk) — )■ H{Ajk{z)) is uniformly bounded, 

fl5.26p implies the norm estimate (see Lemma 15. 4p 



For this estimate to be small when /i — !■ 0, we require the ratio Mq/M to be small enough 
to ensure the condition 

^0 ^^max > ^o/2 > . 

(The bounds fl3.3ip and Mq < Mi show that the ratio Mq/M can indeed be chosen arbitrary 
small.) 

On the other hand, xj -^jki^) is microlocalized inside neigh ( S'j, S'j), so Lemma 15^61 
implies that (Ilj — l)xj -Mjkiz) v'^^ = OHg.{h°°). Putting these estimates altogether, we 
find that 



l^+l 



(5.28) \/v'l G -Hk, M.kiz) v1 = n, M,k{z) v1 + 0(/i^"°/') 

This crucial estimate shows that the projection of Aijk{z) v'l on Hj has a negligible effect. 

We now define the finite rank operators 

(5.29) 

~ dcf jUjMjkiz) Uk : Uk Uj , J G J+{k) , . . . 
Mjk{z) = <■'■' ^ ' m short M(2;) = n/iAl(z) Ilh . 

otherwise , 
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This operators satisfy the same norm bounds (14.361) as their infinite rank counterparts. 
Using these operators, and remembering that the operators R''_ : — ?■ Hg{X) are bounded 
by C(/i-*^o'), we rewrite (K25\i as 

jeJ+(k) 

Generahzing the initial data to arbitrary v+ G "Hi x ■ ■ ■ x "Hatj we obtain the 

Proposition 5.7. Assume z G TZ{6, Mq, h) . Let Vj^ G "H. Then there exists {u,U-) G 
H^{X) X n such that 

(5.30) {i/h){P-z)u + R4z)u^ = 0{h^''>/^)\\v+\\n in Hg{X), 

(5.31) R+{z)u = v+ + 0{h°°)\\v+\\n mK, 

(5.32) MlHaix) < /.-^^°(*— +^) < h~'''^'--\Mn. 
The second part of the solution, is of the form 

M„ = {M{z) - Id)v+ , \\M{z)\\n-,H < h-''°'-- , 
where M{z) = {Mjk{z))j^k=i,..,N is the matrix of operators defined in (I5.29p . 

We collect some properties of the operators Mjk{z), j G J+{k), for z G [—5,5]: 

• Mjk{z) is uniformly bounded, and WF'fj^{Mjk{z)) C Sj x Sk- 

• take pk G Sk, Pj = Fjk,z{pk) e Sf 

(1) if the trajectory segment connecting the points K,k,z{pk), i^j,z{pj) is contained 
in W, then microlocally near {pj,pk), Mjk{z) is an /i-Fourier integral operator 
of order zero with associated canonical transformation Fj^^z = f^jl ° Fj^^z ° i^k.z 

(2) if furthermore the above segment is disjoint from the support of G, then Mjk{z) 
is microlocally unitary near {pj,pk). 

(3) if, on the opposite, this segment contains a part outside W, then there exist 
Xj e C^(neigh(pj)), Xk £ C^(neigh(pfc)), equal to 1 near pj and pk respectively, 
and a time t(pfc) > independent of the exponent M, such that 

XjM,k{z)x';: = 0{h'''^P>=)) : H,^ ^ H,^ . 

For z G 7^(5, Mo, /i) similar statements hold, modulo the fact that the symbol of the Fourier 
integral operator is multiplied by exp{—izt_^./h), which modifies the order of the operator. 

5.3.2. A well-posed inhomogeneous problem. Let us now consider the inhomogeneous prob- 
lem 

(5.33) {i/h){Pg^R- z)u + R^{z)u^ = v v e Hg{X) . 
We will use a partition of unity to decompose v into several component. 
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Take ips E S{T*X), ips = 1 near p~'^ {[-6/2, 6/2]), and ips = outside p-\[-6,6]). The 
operator {Pe^R — z) is elliptic outside p^^[—6/2,6/2]. Taking ips similar with %ps but with 
supp?/'^ C (5/2, 5/2]), the operator 

L iPe,n - z - zz^r) ■ Hi ^ Hg 
is invertible, with uniformly bounded inverse L^^ G \E'^. Hence, by taking 

u = {h/t)L~^{l-i;^)v, 

we find 

{^/h){Pe,R - z)u = {i/h){Pg,R - z - %^^,)u + \\u\\ = (1 - v + 0{h'^)\\v\\ , 

which solves our problem for the data (1 — ip^) v. The first equality uses pseudodifferential 
calculus and the fact that ip^ = 1 on the support of ips: 

Let us now consider the data (ipfv) microlocalized in p~^{[—6,6]). We split this state 
using a spatial cutoff tpR G C^{X), such that tpR = 1 in -6(0, R), ipR = outside -8(0, 2R). 
To solve the equation 

(5.34) {t/h){Pe^R-z)u = v, v = (1 - iPr) iPf v , 
we take the Ansatz 

(5.35) u = E{z)v, 

with E{z) the parametrix of f l4.15p (with P replaced by Pe,R), for the same time T = tmax+e 
as in f l4.39p . It satisfies 

(5.36) (i/h) {Pe^R ~z)u = v- e-*^(^».«-^)/'^ v . 
The time T is small enough, so that 

¥ {p-\[-6, 6]) \ T*(o,«)X) n T^(o,K/2)X = , < t < T . 
Hence, the states 

v{t) = e-^*(^«'«-^)/''^y 

are all microlocalized outside T^(^o r/2)^ ^ ^ [O?^]- The estimate f l3.30p (adapted to the 
weight Go) then implies that [311 Lemma 6.4] 

dtimml^ = ^lm{{Pe,R - z)v{t),v{t))Ha < {-M,/Ci + 2Mo) log(l//i) , Vt G [0,T] , 

where Ci > is independent of the choice of Mi. Once more, we assume Mq/Mi is small 
enough so that -Mi/Ci + 2Mq < -Mi/2Ci, and hence 

SO the problem fl5.34p is solved modulo a remainder 0{h'^^^'^^'^^^). 
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We now consider the component (ipRipfv) microlocahzed in T^^^^^r) ^ P~^{[~^^^])- We 
spht it again using a cutoff ipvi ^ C^(Vi), ipvi = 1 in the set V( d Vi (see the discussion 
after Lemma [5.3p . To solve the problem for the inhomogeneous data 

we use the Ansatz f l5.35p . resulting in the estimate (I5.36p . The microlocalization of u outside 
of V(, together with the assumption (15. lip , implies the norm estimate (see Lemma [5. 4p 



|g-iT(Pg,«-^)/h~ 



Again, we assume Mq/M small enough, so that Mti/2 — MqT > Mti/3. We have solved 
the problem for w up to a remainder 0{h'^^^^^)\\v\\HG- 

We finally consider the data v = ipy-^ipRip'^v microlocahzed inside Vi. For this data, we 
can use the micro local analysis of §4.2.31 If WF/j('i;) is contained inside Vi fl TJ^T ■, then 

WF/i(x^ i?(z) ?;) (see the Ansatz f l4.40p ) will intersect Sj inside the arrival set Ajk{z), but 
not necessarily inside Sj. However, the same phenomenon as in Lemma [5.51 occurs: there 
exists a time ts > such that, for any z G [—5, 5] and any p G Vi fl T~fr{z), 

(5.37) p4p) G S,(^) \ K,,,iS'^) =^ Go(p+(p)) - Go(p) > h ■ 

If we decompose R\ „^{z)E{z)v using the cutoff Xj of (I5.27p . the property f l5.37p implies 
that 

11(1 - xJ) RU{z)E{z)v\\h,^ = 0{h'''^/'-''''^)\\v\\Ha ■ 
Again we assume Mq/M small enough, so that — MqT > Mt^/3. Hence, if we set 



u 



Riiz)xJE{z)v 



= Ri^^{z)xJE{z)v + 0{h^)\\v\\Ha 
= i?i,„(^)E(z)^y+0(/.^^^*3/3)||~||^^^ 

we end up with a solution of fl5.33p modulo a remainder 0{h^^^^^)\\v\\HG- 

We recall that Mi/M is bounded by (13.311) . so all the above error estimates can be put 
in the form (9(/i'^*^)||i;||/^g, with c > independent of M: we have thus shown that the 
problem (15.331) admits a solution for any v G Hq, up to this remainder. We may then 
apply Proposition 15 . 71 to solve the resulting homogeneous problem, and get an approximate 
solution for the full problem (15.221) . We summarize this solution in the following 

Proposition 5.8. Assume z G 7l{6,Mo,h). Let {v,v+) G Hq x Ti. Then there exists 
{u,U-) G Hq X "H such that 



(5.3 



{i/h){P - z)u + R.{z)u. =v + 0{h''^'){\\v\\Ha + \\vA\h) in Hg{X) 
(5.39) \\u\\^. + ll^^ll^ < /,-M.(2W+2e) (ii^ii^^ + ii^^ii^) _ 



50 S. NONNENMACHER, J. SJOSTRAND, AND M. ZWORSKI 

5.4. Invertibility of the Grushin problem. We can transform this approximate solution 
into an exact one. The system fl5.38p can be expressed as an approximate inverse of V{z): 

(5.40) \u~J \v+ 

V{z) £{z) = 1 + n{h) : HgX-H^ HgX-H, \\n{h)\\ = Oih"^) . 

For h small enough the operator I + TZ[h) can be inverted by a Neumann series, so we 
obtain an exact right inverse of V{z), 

£{z)=£{z) (/ + 7^(^))"^ 

Since V{z) is of index zero, S{z) is also a left inverse, which proves the well-posedness of 
our Grushin problem f l5.22p . 

Theorem 2. We consider h > small enough, and z G TZ{6, Mq, h) . For every {v,v^) G 
Hg X Ti, there exists a unique {u,uJ) G Hq x "H such that 



(5.41) 



(i/h){Pe^R — z)u + R-{z)u- =v 
Rjf-{z)u = f+ 



where R±{z) are defined by (I5.18P and fl5.20p . The estimates liS. 39\) hold, so if we write 

£(z)-- 





1 = 




: 











E 

E_ E_^ 

then the following operator norms (between the appropriate Hilbert spaces ) are bounded by: 

(5.42) \\E\\, \\E+\\, \\E^\\, = 0(/i-*^«(2*--+2^)). 
Moreover, we have a precise expression for the effective Hamiltonian: 

(5.43) E.^{z) = -I + M{z) + Ou^nih^''') = -I + M(z, h) , 

where M[z) is the matrix of "open quantum maps" defined in (15.290 and described after 
Proposition \5. 7\ 

Remark. If we restrict the parameter 2; to a rectangle of height | lmz\ < Ch instead of 
llm^l < Mo/ilog(l//i), the bounds (ICTD become \\E4z)\\ = C(l). 

Theorem [Hand the formula (11.80 follow from this more precise result. In fact, the equality 
(I3.13P shows that 

(5.44) rank (p xR{w)x dw = rank ® x^d r{'^)x dw = tr ® Rg r{w) dw , 

Jz Jz ' 27rz ' 

see [m Proposition 3.6] for the proof of the last identity in the simpler case of compactly 
supported perturbations, and [39], Section 5] for the general case. 
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The well-posedness of our Grushin problem means that we can apply formula fl3.33p 
recalled in §3.71 It shows that the right hand side in fl5.44p is equal to 

— tr (p E^+{w)-^E'_Jw) dw, 

which in view of (15.43^ gives ( ll.Sp . The exponent L =^ c'M in the remainder of ( I5.43P 
depends on the integer M > used in the scaling of the weight function G, which can be 
chosen arbitrary large, independently of c' > 0. 
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